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ABSTRACT 

We have obtained high-dispersion spectra for a sample of 256 pre-main-sequence 
stars in the Orion Nebula Cluster (ONC) in order to measure their projected rotational 
velocities and investigate the rotational evolution and physical properties of young, low- 
mass stars. Half the stars were chosen because they had known photometric periods and 
the other half were selected as a control sample of objects without known periods from 
the same portion of the H-R diagram. More than 90% of the spectra yielded v sin{i) 
measurements, although about one-third are upper limits. We find strong evidence 
confirming the long-held assumption that the periodic light variations of T Tauri stars 
are caused by rotation of spots on their surfaces. We find no statistically significant 
difference between the v sin{i) distributions of the periodic and control samples, indi- 
cating that there is no strong bias in studying the rotational properties of young stars 
using periodic variables. Likewise, the classical and weak T Tauri stars exhibit v sin{i) 
distributions that are statistically the same. For stars with known period and vsin{i), 
the mean value of sin{i) is significantly lower than expected for a random distribution 
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of stellar rotation axes. This could be caused by errors in one or more of the quantities 
that contribute to the sin{i) calculation or to a real physical effect. Wc investigate the 
possible causes and find that {sin{i)) has the expected value if wc increase the effective 
temperatures of our stars by 400—600 K. Finally, we have calculated minimum radii 
{R sin{i)) for stars with both v sin{i) and period, as well as average radii for objects 
grouped by their location in the H-R diagram. We find evidence at the three-sigma 
level that the radii of the stars on similar mass tracks are decreasing as the stars move 
closer to the ZAMS. 

Subject headings: stars: pre-main sequence; stars: rotation; stars: fundamental param- 
eters 

1. Introduction 

T Tauri stars are low-mass, late-type stars that are still undergoing gravitational contraction 
toward the main sequence. They often show characteristics that reflect the dynamic, turbulent 
nature of this early stage of evolution: powerful magnetic fields cause prominent dark spots to 
form on their surfaces, bipolar outflows or jets direct material away from them at high velocities, 
and accreting matter that is channelled along magnetic field lines produces surface hot spots that 
give rise to rapid brightness fluctuations (Bertout 1989; Herbst et al. 1994). In addition, many 
T Tauri stars have thin circumstellar disks of dust and gas that can extend out to a few hundred 
astronomical units. Studying the rotation and evolution of these objects can lead to insights into 
the early evolution of solar-type stars and the planetary systems that form around them. 

One of the best regions of the sky in which to study T Tauri stars is the Orion Nebula Cluster 
(ONC). The ONC is located in the portion of the Orion Nebula centered on four OB-type stars 

called the Trapezium. The cluster is especially well-suited for investigations of star formation and 
early stellar evolution because of its large stellar population, the range of stellar masses present 
(~0.1 - 50 Mq), and the very young ages (< 1-2 Myr) of its member stars (Hillenbrand 1997). 
In addition, cluster membership for nearly 900 stars in the ONC has been established by a proper 
motion survey by Jones Sz Walker (1988), and Hillenbrand (1997) and Hillenbrand et al. (1998) 
have published extensive studies of the observable properties of ONC members. These reasons, 
combined with the relative proximity of the ONC (its distance is ~470 pc; Walker 1969), make it 
a natural choice for detailed studies of the physical and rotational properties of T Tauri stars. 

We began such a study in 1997, with the aim of using high-dispersion spectroscopy to measure 
projected rotational velocities, or v sin{i) values, for a large sample of pre-main-sequence stars 
in the ONC. Rotation periods, measured via starspot-modulated light curves, have already been 
obtained for about 400 stars in the ONC (Herbst et al. 2000; Stassun et al. 1999; Herbst et al. 
2001). We sought to combine data obtained via both photometric and spectroscopic techniques 
in order to investigate the fundamental properties and rotational evolution of young, low-mass 
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stars. When we began our study, periods for approximately 130 ONC stars had been measured at 
Van Vleck Observatory (VVO) on the campus of Weslcyan University (see Mandel & Herbst 1991, 
Attridge & Herbst 1992, Choi &: Herbst 1996). These will hereafter be referred to as the VVO 
periods. Details about the VVO observing program are described in Herbst et al. (2000); briefly, 
multiple fields around the ONC are imaged repeatedly on each of ~30 clear nights over the course 
of an observing season, using a 0.6-meter telescope and a CCD. For our v sin{i) study, we chose 
as targets a sample of stars with measured VVO periods (hereafter the "Periodic sample"), as well 
as a comparably-sized sample of ONC stars from the same portion of the H-R diagram without 
known rotation periods (hereafter the "Control sample"). The properties of these two samples are 
described in more detail in Section 2.1. 

In addition to using our v sin{i) measurements to investigate pre-main-sequence stellar rotation 
in the ONC, we had other specific objectives in mind. First of all, measuring v sin{i) values for stars 
with known periods can confirm that rotation of surface cool or hot spots is indeed the cause of the 
periodic variations of T Tauri stars. It has been assumed since the pioneering studies of Rydgren 
and Vrba in the 1980's (e.g., Rydgren & Vrba 1983, Vrba et al. 1986) that the periodic variations 
of these objects can be interpreted as due to rotation. Stassun ct al. (1999) observed a correlation 
between line width and angular velocity among ONC stars, supporting this assumption. We address 
the issue further in this paper. Also, for individual stars there is often an ambiguity of a factor of 
two in the rotation period (i.e., caused by spots on two sides) as well as the possibility of aliasing. 
Combining v sin{i) and period measures for a sample of stars of known luminosity and effective 
temperature can allow us to test the rotational interpretation of periods for pre-main-sequence stars 
and can help us identify cases of period doubling or aliasing. 

Furthermore, for objects with measured rotation periods, we can use their v sin{i) values to 
obtain statistical estimates of their radii. With the exception of four known pre-main-sequence 
eclipsing binaries (Popper 1987; Casey et al. 1995, 1998; Andersen 1991; Mamajek et al. 2000; 
Covino et al. 2000), we have no direct measurements of the sizes of pre-main-sequence stars. Since 
Pvsin{i) = 2iTRsin{i), combining the stellar rotation period P with v sin{i) can constrain the 
stellar radius R to within a factor of sin{i). The average value of sin{i) can be calculated; thus 
average radii can be derived for stars grouped according to their location in the H-R diagram. In 
combination with luminosities, radii also provide an independent check on effective temperatures, 
which are typically derived from spectral types by comparing pre-main-sequence stellar spectra 
with those of main sequence or evolved stars. 

Another objective of the study was to determine whether the sample of stars that have per- 
sistent spots (and therefore can yield rotation periods) have rotational properties that are repre- 
sentative of the properties of ONC members as a whole. One can imagine that our understanding 
of pre-main-sequence rotational evolution might be biased by data derived from spotted stars. For 
example, if rapidly-rotating stars have more active dynamos, they may be more likely to have 
prominent, persistent spots on their surfaces and thus would be overreprcscnted in photometric ro- 
tation studies. Comparing the v sin{i) distributions of the Periodic and Control samples provides 
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an important observational test of our understanding of the initial conditions of pre-main-sequence 
rotation. 

Here we present the results from our v sin{i) study of 256 targets in the ONC. The following 
section contains a summary of the characteristics of the stellar sample and describes the observations 
and technique used to measure v sin{i). The results of the study are presented in Section 3. In 
that section, we explore the measured v sin{i) values of various sub-samples of stars, compare our 
vsin{i) values to the equatorial velocities calculated for stars with measured period, luminosity 
and effective temperature, and finally combine period and v sin{i) data to derive average radii for 
stars binned according to their location on the H-R diagram. A summary of our findings is given 
in the last section of the paper. 

2. Observations and Analysis 

2.1. Characteristics of the Sample 

For information about our target stars, we have relied on two previous studies of the ONC. 
Jones k, Walker (1988) undertook a proper motion study of the cluster and published membership 
probabilities for 1053 stars in the region. Throughout this paper, we refer to stars by the number 
assigned by Jones & Walker, e.g., JW 100. Hillenbrand (1997) published a comprehensive study 
of star formation in the ONC, including optical photometric data for --^1600 stars and spectra for 
~900 stars. All luminosities, effective temperatures, spectral types, / magnitudes, and V—I colors 
given in this paper come from the Hillenbrand study. 

The initial set of stars observed for this work was selected because the stars had measured 
WO periods. After the first observing run (see Section 2.2), a sample of stars without known 
periods — the Control sample — was chosen. In order to make a useful comparison between the 
V sin{i) distributions of the Periodic and Control sample stars, the latter were chosen to match 
the basic characteristics of the Periodic sample. The properties of the overall stellar sample are as 
follows (the few stars with characteristics outside those of the larger sample are noted): (1) Spectral 
types in the range G6 to M6, with the large majority of stars (88%) in the range KO to M4. There 
are three exceptions: JW 945 (spectral type B6), JW 1033 (spectral type AO— A5), and JW 165 
(spectral type A7). (2) I magnitudes from 10—16, and V—I colors ranging from 0.8 to 4.3. (3) 
Typical masses (as determined by Hillenbrand (1997) from prc-main-scqucncc evolutionary tracks) 
in the range 0.1 to 3 M©; the median mass is 0.3 Mq. (4) Cluster membership probabilities from 
Jones & Walker (1988) equalling 79% or above for all but four stars. The large majority of stars 
(78%) have membership probabilities of 99%. Four stars were observed that have low membership 
probabilities; two of them (JW 308 and JW 794) had been found to be periodic in the WO study 
at the time they were observed, and the other two (JW 3 and JW 1021) have been found to be 
periodic in a subsequent study (Hcrbst et al. 2001). All four stars exhibit the low-amplitude, 
sinusoidal light variations typical of T Tauri stars. It is possible that they are T Tauri stars within 
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the ONC and that their proper motions are anomalous or in error, or that they are T Tauri stars 
superimposed on the ONC but not members of the chister. For three of them (JW 3, JW 308 
and JW 1021) we measure hchocentric radial velocities that are within 1—2 standard deviations 
from the mean of our distribution (at 27 kms~^, 20 kms~^ and 35 kms~^, respectively) indicating 
that these may indeed be ONC members. For JW 794 we measure a radial velocity of 124 kms~^, 
indicating that this may be a true non-member of the cluster. 

Because we wished to have roughly equal numbers of objects in the Periodic and Control 
samples, and WO periods were known for ~130 stars, we needed observations of a total of ap- 
proximately 260 targets. We were able to observe 118 stars for the Periodic sample and 138 for 
the Control sample. Table 1 gives basic information for the 256 stars observed. Columns 1 and 3 
give the identification number and membership probability from Jones & Walker (1988). Column 2 
specifies whether the star was originally chosen as part of the Periodic (P) or Control (C) sample. 
Columns 4, 5, 6 and 7 are, respectively, /, V—I, the log of the luminosity, and the log of the effective 
temperature for the star from Hillenbrand (1997). Columns 8 and 9 are the measured near-infrared 
excess and the equivalent width of the Ca II A8542 line from Hillenbrand ct al. (1998). Column 
10 is the measured period of variation of the star from photometric data, and column 11 indicates 
the reference from which the period was taken (V = VVO data (Herbst et al. 2000), S = data 
from Stassun et al. (1999), E = ESO data from Herbst et al. (2001)). Column 12 gives the star's 
V sin{i) and estimated error measured in this study, and column 13 lists the projected radius of 
the star, calculated by combining the v sin{i) and rotation period (see Section 3.2). Column 14 
indicates whether the cross-correlation function showed sufficient structure to indicate that the star 
may be a double-lined spectroscopic binary (see Sections 2.4 and 3.1 for details). The last column 
in the table specifies when the star was observed. Stars that initially showed structure in their 
cross-correlation functions were, when possible, observed more than once. 

The relative spatial locations of the stars observed for this survey compared to the overall Jones 
& Walker (1988) survey are shown in Figure 1. Small dots in the figure show the locations of the 
1053 stars surveyed by Jones &; Walker, and filled circles are the 256 stars observed for this study. 
The observed sample is spread throughout the area surveyed by Jones Sz Walker. The upper left 
portion of the field that is relatively empty of stars is a region of the ONC that is optically obscured. 
Figure 2 shows the positions of the target stars on an HR diagram. Stars in the Periodic sample 
are shown as filled circles and Control sample stars appear as open circles. The pre-main-sequence 
evolutionary tracks shown in the figure come from D'Antona & Mazzitelli (1994, model 1). The 
target stars span the mass range of the tracks and most lie above the 3xl0^-year isochrone. The 
bulk of the stars are concentrated toward the low-temperature, low-luminosity end of the diagram. 

We note that the Control sample stars were chosen without regard to whether they appeared 
within the dozen or so fields that had been monitored in the WO program. This means that they 
are not necessarily non-periodic stars (i.e., we were not selecting against stars with measurable 
periods), but are simply ONC members with the same basic characteristics as the Periodic sample 
stars. In that sense they were selected to be a true "control" sample. Since this study was initiated. 
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about 30% of the Control sample stars have had their periods measured. Thus Table 1 contains 
stars that are marked with a "C" in the Sample column, but also have a measured period listed. 

2.2. Observations 

Because we needed to obtain high-dispersion spectra for a large number of objects over a fairly 
small region of sky (the ^1/2 degree subtended by the ONC) a multi-object spectrograph seemed 
the ideal instrument for our purposes. The 3.5-meter WIYN telescope^ at Kitt Peak National 
Observatory is equipped with the Hydra Multi-Object Spectrograph, which makes possible simul- 
taneous spectroscopy of ~100 targets over a 1-degree field. Observations for this study took place 
at WIYN during three separate observing runs, in January 1997, December 1997 and December 
1998. 

For all observations we used the Bench Camera, Red Fiber cable, and 316@63.4 Echelle grating. 
This spectrograph setup gave us fairly high resolution {R ^ 21,500) while still allowing us to observe 
our faintest targets {V ~ 18) with reasonable integration times. The spectra cover the region 6240 
A to 6540 A, with Xcentral = 6400 A, and a dispersion of 0.145 A per pixel. We chose to work in 
the red portion of the spectrum because of the colors of our program stars {V—I = 1—4), while at 
the same time avoiding the region of the spectrum around Ha, where nebular emission can easily 
overwhelm the light from the star. The spectral resolution (corresponding to a typical FWHM of 
the slit profile of 2 pixels) is 13.6 kms~^ at 6400 A. 

The 256 target stars were observed in five different fiber configurations, so typical configurations 
had fibers positioned on ~50— 60 objects. Crowding in the ONC and limitations on exactly how 
Hydra fibers can be positioned prevented us from observing larger numbers of stars in a given 
configuration. At least 15—20 of the remaining fibers in each pointing were placed on blank sky, 
to be used for sky subtraction during processing of the spectra. Total integration times for the 
program star fields ranged from ~3— 4 hours, depending on the magnitudes of the stars in the 
configuration and the sky conditions. In addition to the program stars, we observed five objects 
from the Gliese Catalog of Nearby Stars (Gliese & Jahreiss 1991) with spectral types spanning the 
range of our targets; the spectra of these stars were used as narrow-lined templates in the cross- 
correlation process. Information about the templates is given in Table 2. Daytime sky spectra were 
also obtained for use as G2V template spectra. 

Figure 3 shows examples of narrow- lined template and object spectra obtained for this study. 
The spectra are arranged in order of increasing rotational velocity, and illustrate the effect that 
this has on stellar absorption lines. The top spectrum is that of Gliese 144, a K2 dwarf used 
as a narrow-lined template. Its rotational velocity is below our v sin{i) measurement limit, and 
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its lines appear narrow. The spectra in the bottom four panels are target stars with projected 

velocities ranging from ~20— 100 kms~^. The spectral lines become broader and shallower with 
larger v sin{i). The emission lines present in a few of the spectra are nebular lines at A6363.8 due 
to [O I] and A6371.3 due to Si II. Regions of the spectra with emission lines were excluded from the 
vsin{i) measurement process (i.e., not included in the cross-correlation) described in Section 2.4. 

2.3. Initial Reductions 

Preliminary data reduction was accomplished with standard IRAF^ tasks, e.g., ZEROCOM- 
BINE to create combined bias images, CCDPROC to bias-subtract and trim the object images, 
dome flats and comparison lamp frames, and FLATCOMBINE to create combined dome flats. Then 
the IRAF task DOHYDRA was used to perform spectral extraction, flat-fielding, fiber throughput 
correction, wavelength calibration, and sky subtraction on the object frames. The wavelength cali- 
bration was done using observations of a Thorium- Argon comparison lamp. In the sky subtraction 
step, sky fibers with unusually high signal (due to contamination by an object) were identified 
and deleted. A combined sky spectrum was then made by averaging the remaining sky spectra 
and using a sigma-clipping algorithm to eliminate cosmic rays. We note that the sky fibers placed 
randomly around the field will collect flux from both the sky and the Orion nebula. The sky values 
varied spatially, reflecting variations in the nebular continuum. Although the sky level in the vast 
majority of fibers varied by 1 10%, a few had deviations of as much as 30%. A series of tests 
was performed that indicate that a variable background does not affect the final measured v sin{i) 
values. The template star and daytime sky spectra were reduced in the same manner as the target 
star spectra, with the only difference being that sky subtraction was not executed on the daytime 
sky frames. 

Once DOHYDRA had been run on the object frames, individual integrations taken of the same 
field (within a given observing run) were scaled and then combined into a single frame using the 
task SCOMBINE, to remove cosmic rays and produce a single, high signal-to-noise observation. 
Because the signal degrades at the ends of some of the images, the frames were clipped so that the 
spectra cover the region between 6275 and 6525 A. All stellar spectra were continuum-normalized 
before the cross-correlation step was executed. 

2.4. Measuring Radial and Rotational Velocities 

Rotational and radial velocities were derived using the IRAF task FXCOR, which performs 
Fourier cross-correlation of an object spectrum against a specified template. Figure 3 shows that 
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increasing rotational velocity causes stellar absorption lines to become broader and more shallow. 

When a stellar spectrum with rotationally-broadcned lines is cross-correlated against a narrow-lined 
spectrum, the width of the cross-correlation peak is sensitive to the amount of broadening present 
in the first spectrum. Thus, by measuring the width of the peak, one can obtain a measurement of 
the rotational velocity of the star. 

To measure v sin{i) for our target stars, we began by calibrating the relationship between the 
width of the cross-correlation peak and v sin{i). To do this, a narrow-lined template spectrum was 
artificially "spun up" to mimic that of a higher-velocity star by convolving it with a theoretical 
rotation profile (Gray 1992). A series of such broadened spectra was created for v sin{i) values 
spanning the expected range of observed values. Each broadened spectrum was then cross-correlated 
against the original narrow-lined spectrum, and the FWHM of the cross-correlation peak was 
measured. The resultant relationship between v sin{i) and the FWHM of the cross-correlation 
peak was fit with a 4th-order polynomial. This process was executed for each of the six template 
stars (including the Sun). Figure 4 shows, as an example, the calibration curve data produced by 
cross-correlating broadened versions of the spectrum of the template star Gliese 114 against the 
unbroadcncd spectrum. The points on the figure are the data and the dotted line is the polynomial 
fit to the data. Below 7 kms^^, the relationship between the FWHM of the cross-correlation peak 
and V sin{i) is flat, but from 7 kms^^ upward it is fairly linear. Figure 4 illustrates the absolute 
lower limit on a sin{i) measurement using our analysis method, under ideal circumstances. 

The actual limit to which we can measure v sin[i) depends in part on the size of the slit (fiber) 
image. With the Bench Camera and our spectrograph setup, the slit image width (as determined 
by measuring the FWHM of emission line profiles in the comparison lamp spectra) varied from, 
typically, 1.5 to 1.7 pixels toward the ends of the spectra to >2 pixels in the middle (in the worst 
cases, small portions of some of the spectra have FWHM as large as 2.8—2.9 pixels). This makes 
determining the exact limit on v sin{i) measurement more difficult than if the focus were uniform 
across the CCD chip. In the regime in which the broadening of a spectral line can be measured 
(i.e., once the width due to intrinsic broadening is larger than the slit width), v sin{i) is related 
to the half-width at zero intensity (HWZI) of the line. The base of the line, or full-width at zero 
intensity (FWZI), is on the order of two times the FWHM of the line, which means that HWZI 
is approximately equal to FWHM. The FWHM of the slit image varies in our spectra from <2 
pixels to >2 pixels. Spectral lines from all of these regions contribute to the cross-correlation. In 
places where, for example, the FWHM is 1.7 pixels, the base of a spectral line spans 3—4 pixels 
and therefore FWZI is resolved. In such circumstances v sin{i) values somewhat below the 2-pixel- 
FWHM spectral resolution of 13.6 kms~^ can be measured. 

Because of the focus variation we cannot determine with certainty what the limit on measuring 
V sin{i) is on a star-by-star basis. Our calibration curve suggests that we are sensitive to rotational 
broadening as small as 7 km s~^. However, without observations of v sin{i) standards in that regime, 
we are not comfortable reporting v sin{i) measurements much below our spectral resolution. Thus 
we take our limit on v sin{i) measurement to be 11 kms~^, since this corresponds to a best-focus 
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value of between 1.6 and 1.7 pixels (FWHM). We emphasize that reported values at or somewhat 
above this limit should be treated with caution. In some cases these values could be upper limits 
rather than actual velocity measurements. 

Once the relationship between the width of the cross-correlation peak and v sin{i) was es- 
tablished for the narrow-lined templates, each program star spectrum was cross-correlated against 
one of the template spectra. The template star closest in T^s to the object was used as the cross- 
correlation template. The FWHM of the cross-correlation peak was measured and the v sin{i) value 
for the target star calculated using the derived relationship between those two quantities. Radial 
velocities for the target stars were obtained as part of the cross-correlation step. Radial velocities 
used for the template stars come from the Gliese catalog (Gliese & Jahreiss 1991). 

In addition to being useful for measuring radial and rotational velocities, the cross-correlation 
function can also show structure that may indicate that a star is a double-lined spectroscopic 
binary (SB2). The peak of the function, instead of appearing Gaussian in shape, may look like 
two Gaussians superimposed on each other, due to the presence of two sets of lines in the stellar 
spectrum. Figure 5 shows example cross-correlation functions for two stars. The top panel of the 
figure is the peak of the cross-correlation function for JW 790. The peak is fairly regular and 
nearly Gaussian in shape. The bottom two panels are the cross-correlation function peaks for 
JW 669, made using spectra from January 1997 and December 1997. The structure in this star's 
cross-correlation peak indicates that it may be an SB2. 

3. Results 

3.1. Rotational Velocities 

Of the 256 stars observed for this study, 18 did not yield reliable v sin{i) values, primarily 
because their spectra lacked sufficient signal. In these cases, there was either no obvious peak 
present in the cross-correlation function, or the peak found by the FXCOR task was not far enough 
above the noise to be believable. Eighty-three of the remaining 238 stars have v sin{i) less than 
or equal to our estimated limit of 11 kms~^. The measured v sin{i) values for our program stars 

are given in Table 1. For the 18 stars for which no v sin{i) was measured, " " is given in the 

V sin(i) column; for stars with v sin{i) below our estimated limit, "<11.0" is given in that column. 

Seven stars had sufficient structure in their cross-correlation functions to indicate that they 
may be SB2s. The cross-correlation peaks for these stars are shown in Figures 5 and 6. Choosing 
SB2 candidates was a fairly subjective process; we simply examined the cross-correlation peaks and 
noted those that were shaped more like two (or more) Gaussians superimposed than like a single, 
smooth Gaussian curve. Other stars also showed structure in their cross-correlation peaks, in the 
sense that the peaks appeared noisy and irregular rather than smooth, but we did not mark a star 
as an SB2 candidate unless its peak looked like multiple Gaussians. The seven SB2 candidates 
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are indicated in Table 1 with a "Y" in the "SB2?" column. If any of them are in fact SB2s, their 
V sin{i) values would likely be too large, since some of the width of the cross-correlation peak would 
be due to the fact that their spectra contain two sets of lines. 

We have calculated uncertainties on our v sin{i) values by making use of the parameter r from 
Tonry & Davis (1979). It provides a measure of the signal-to- noise of the cross-correlation peak; 
specifically, r is the ratio of the height of the cross-correlation peak to that of the average peak 
in the noise component of the cross-correlation function. Tonry &; Davis (1979) showed that the 
errors in velocities determined via cross-correlation should be proportional to (1 -|- r)~^. Hartmann 
et al. (1986) used cross-correlation to measure v sin{i) for stars in Taurus- Auriga and Orion, and 
investigated their observed errors as a function of r. They found that the quantity ±v sin{i)/(l + r) 
provided a good estimate for the 90% confidence level of one of their v sin{i) measurements. 

We have adopted the value zizvsin{i)/{l +r) as a reasonable esfAm,ate of the onc-sigma errors 
on our V sin{i) values. We can check this by comparing v sin{i) values for the eight stars for which 
we have repeated observations. Two of the stars, JW 2 and JW 433, have v sin{i) values below 
our estimated measurement limit. For the remaining six stars, the percentage difference between 
the V sin{i) measurements from different observing runs ranges from 4% to 24%, with an average 
difference of 13%. (We note that two of these six stars, JW 99 and JW 669, are SB2 candidates, so 
the widths of their cross-correlation peaks may change due to orbital motion of the member stars, 
if they are indeed binaries. However, the percentage differences between the v sin{i) values from 
their observations, at 14% and 4%, respectively, are consistent with the differences for the other 
stars.) The relative v sin{i) errors calculated using the r parameter for the six stars range from 10% 
to 26%, with an average error of 20%. Furthermore, when we compare the observed percentage 
difference in v sin(i) to the relative v sin{i) errors calculated from r on a star-by-star basis, we find 
that the uncertainty calculated from r is equal to or larger than the observed difference for all six 
of the stars. Thus the average error calculated from r is comparable to or larger than our estimated 
precision based on repeated measurements. The quantity ±vsin(i)/(l + r) is given as the error 
on V sin{i) in Table 1. For the stars with repeated measurements, the v sin{i) and error given in 
Table 1 are the weighted mean and its associated error, calculated from the multiple measurements. 

In the absence of a larger number of cases in which stars had repeated observations taken 
during separate observing runs, we can perform another test of whether our v sin{i) error estimates 
are appropriate. As explained in Section 2.3, multiple observations of a given Hydra field taken 
during the same observing run were combined to remove cosmic rays and create a single, high signal- 
to-noise observation. A useful comparison to make is to take these individual Hydra frames and 
measure v sin(i) for the stars appearing in them. This is equivalent to having multiple, independent 
observations of the stars in these fields. Moreover, the signal-to-noise ratios of the spectra in each 
frame will differ depending on the observing conditions during the individual integrations, providing 
an additional test of how well our v sin{i) measurements and errors hold up under varying levels of 
signal-to-noisc. To perform this test, we analyzed Hydra frames taken of two different fields, one 
observed during December 1997 and another in December 1998. In the former case, six individual 
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integrations of the same field were taken over a period of about three weeks. In the latter case, 
five exposures of another field were obtained on a single night. To measure v sin{i) in the single 
frames, we used the identical method described in Section 2.4 used to measure velocities in the 
combined frames. Because the integration times for the single frames were only 1/5 to 1/6 of the 
total combined integration, many of the stars in the single frames had insufficient signal to measure 
vsin{i). Sixty-one objects did yield reliable multiple measurements, and of these, 34 have v sin{i) 
> 11 kms~"^. For these 34 stars, we compared their calculated v sin{i) errors {^vsin{i)/{l + r)) 
with the standard deviation of the individual v sin{i) values from the multiple measurements. For 
31 of 34 stars (91%), the v sin{i) error was equal to or larger than the standard deviation of the 
individual measurements, and in a handful of cases much larger. For the remaining three stars, the 
computed v sin{i) error equalled ~60— 80% of the standard deviation. Overall the test indicates 
again that our adoption of rk.v sin{i) /(I + r) as a reasonable estimate of the one-sigma error on the 
V sin{i) values is both valid and appropriate. 

A possible source of systematic error in the v sin{i) measurement process comes from template 
mismatch, i.e., an error that may arise if an object spectrum is cross-correlated against a template 
star that does not exactly match the object's spectral type. Template mismatch might be a factor 
in our data set because we do not have a template star to match each individual spectral type of 
our objects, but instead have a set of five templates spanning the approximate range of our targets. 
To investigate the effect of template mismatch, a set of 40 stars was randomly chosen from the 238 
for which v sin{i) was successfully measured. These 40 stars were cross-correlated against template 
stars hotter than and cooler than the correct template based on the star's spectral type. For 
example, a K6 target that was originally cross-correlated against the K7V template star Gliese 114 
was for this test cross-correlated against Gliese 144 (K2V) and Gliese 15A (M1.5V). Comparisons 
were made of the v sin{i) values resulting from these tests for the 27 stars with v sin{i) values of at 
least 11 kms"-*^. In 22 of the 27 cases (82%), the v sin{i) values measured using the hotter or cooler 
templates matched the original v sin{i) within the one-sigma error bars, and in all 27 cases they 
matched within less than three sigma. The mean difference between the original v sin{i) value and 
that measured using a hotter template was —2.5 kms^^, with a standard deviation of 3.9 kms^^. 
The mean difference between the original value and the value when a cooler template was used 
was -1-1.3 kms~^, with standard deviation 5.1 kms~^. In other words, using a template that is 
too hot for a given target star causes the v sin{i) measured to be too large, and the opposite is 
true for a template that is too cool. In the case of the too-cool template, the mean difference is 
not statistically significant, whereas it is marginally so for the templates that were too hot. We 
note however that the test performed mimics an extreme situation in which the stars have been 
cross-correlated against templates that are several spectral subclasses away from the appropriate 
spectral type. Even under such circumstances, the mean systematic difference for the group of 27 
stars is very small. The overall agreement within the errors for individual stars indicates that the 
uncertainty associated with template mismatch is accounted for by our estimated errors on v sin{i). 

Another potential source of error in our v sin{i) values is line blending in the stellar spectra. 
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Including regions of the spectra that have strong hne blending could bias our results toward larger 
derived vsin{i). To test for this effect, 40 stars were randomly chosen (a different set than used 
for the template mismatch test) from the 238 with measured v sin{i). Several regions with possible 
line blending were identified using a high signal-to-noise daytime sky spectrum. These regions were 
then excluded from the cross-correlation process. The resultant v sin{i) values for 28 stars for 
which V sin{i) was at least 11 kms~^ were compared with those originally measured. For 27 of the 
28 stars (96%), the v sin{i) values from the test were the same as the original values within two 
sigma, and in all 28 cases they agreed within less than three sigma. The mean difference between 
the original v sin{i) and that measured with line-blended regions excluded was 0.45 kms^^, with a 
standard deviation of 3.0 kms^^. This difference is not significant. Thus, as in the case of template 
mismatch, we again find that the potential systematic uncertainties associated with line blending 
are taken into account by the uncertainties we have estimated for our v sin(i) values. 

Figure 7 shows a histogram of the heliocentric radial velocity values for the stars for which 
vsin{i) was measured. The outlier appearing at —8 kms^^ is JW 50, one of our SB2 candidates; 
to derive the radial velocity for this star, we (arbitrarily) chose to fit the blue peak in the cross- 
correlation function, since the two peaks were well-separated. The star with a radial velocity of 
87 kms^^ is JW 363, which has an ONC membership probability from Jones & Walker (1988) of 
98%, so this star may be an SBl. The mean of the distribution, excluding JW 363 and all SB2 
candidates, is 26.7 kms~^, with a standard deviation of 5.6 kms~^. This is comparable (although 
with larger errors) to the expected mean radial velocity for stars in this region, which is ~25 kms~^ 
with a dispersion of ~2 kms~^ (Jones & Walker 1988; R. Mathieu, unpublished data). 

Fourteen of the stars in our sample with v sin{i) >11 kms~^ have had their v sin(i) values 
measured in other studies. Figure 8 shows a comparison between our measured v sin{i) values and 
those from work by Duncan (1993) and Wolff, Strom, &: Hillenbrand (2001). In the figure, squares 
represent stars that our study has in common with Duncan's study, and triangles are stars we have 
in common with Wolff, Strom, &: Hillenbrand. Duncan estimated that his v sin(i) values have an 
accuracy of 10—15%, so the squares are shown with horizontal error bars equal to 15% of vsin{i). 
The points in the figure scatter around the line marking equal v sin{i) values from our study and 
the others; our measured values arc larger for nine of the stars, and smaller for the remaining five. 
The mean difference between the v sin{i) values is 0.5ib2.0 kms^"*^, indicating reasonable agreement 
between our values and those from previous work. The figure does show, however, that the random 
errors of all three studies have probably been underestimated, since about half of the points do 
not intersect the equality line within the error bars. Nevertheless, there are no obvious systematic 
differences between our v sin{i) values and those measured by other groups. 

3.1.1. V sin{i) Distributions 

The uppermost plot in Figure 9 shows the v sin(i) measurement results from this study, in the 
form of a histogram of the v sin{i) distribution for all 238 stars for which a v sin(i) was measured 
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(including upper limits). In this and the other two plots shown in Figure 9, the number of stars in 
each bin has been normalized by the total number of stars in the particular sample. Approximately 
68% of the stars have v sinii) less than 20 kms^-*^. The remaining 32% of v sin(i) values are 
distributed in a "tail" that extends from 20 knis^-*^ to '^lOO kms~^. 

One of the first systematic studies of T Tauri star rotational velocities was done by Vogel &: 
Kuhi (1981), who observed 50 stars in the Taurus- Auriga and Orion complexes. The Vogel & Kuhi 
sample was somewhat brighter than our sample, having V magnitudes in the range 11—15, whereas 
our stars have V ~ 12—19. Their stars were generally more massive than ours, with masses larger 
than ^0.5 M©. Vogel & Kuhi found that three-quarters of the stars they observed had velocities 
below their observational limit of 25—35 kms~^, and concluded that stars with masses < 1.5 
generally have v sin{i) less than 25 kms~^. A 1986 study by Hartmann et al. of 50 stars in Taurus- 
Auriga and Orion found similar results. They observed stars with V magnitudes brighter than 
13—14, thus comparable to the brightest stars in our sample. Hartmann etal. concluded that 
nearly all the stars in the Taurus- Auriga complex have v sin{i) < 20km s~^. Although the large 
majority of our sample stars have v sin{i) < 20 kms~^, 25 stars have v sin{i) > 45 kms~^. Twenty 
of the 25 stars with large v sin{i) values have masses less than 1.2-Mq, and the star with the largest 
V sin{i) value is JW 526, a O.47-M0 star. Thus the finding of Vogel & Kuhi that stars with mass 
< 1.5 M0 have velocities less than 25 kms~^ is not true for our sample. This is consistent with 
photometric studies (e.g., Choi & Herbst 1996, Stassun etal. 1999, Herbst etal. 2000, 2001) which 
show that a number of the low- mass stars in the ONC have short rotation periods (e.g., less than 
a day or two). 

One of the objectives of this study was to test whether the rotation data gained from photomet- 
ric studies of pre-main-sequence stars (i.e., studies in which periods are measured via photometric 

monitoring) are biased in some way, or are truly representative of the rotational properties of 
young stars. To address this issue, wc compared the v sin{i) distributions of the Periodic sample 
(the sample of stars known to be periodic from the VVO monitoring study) and the Control sample 
(the sample of stars selected to have luminosities and effective temperatures similar to the Periodic 
sample stars). (As noted in Section 2.1, about 30% of the Control sample stars have since been 
found to be periodic.) This comparison is shown in the two lower plots in Figure 9, which show 
the V sin{i) distributions for 111 stars in the Periodic sample and 127 stars in the Control sample. 
Both distributions look similar to that of the full sample. One might expect that studies based on 
rotation periods could be skewed toward including a higher fraction of rapid rotators than appears 
in the overall population of young stars, because rapidly-rotating stars with more active dynamos 
might tend to have more prominent and/or persistent surface spots. In that case, the Periodic sam- 
ple should have a larger proportion of rapid rotators compared to the Control sample. This may at 
first glance appear to be true, since there is a slight overabundance of rapid rotators in the Periodic 
sample distribution. (23% of stars in the Periodic sample have v sin(i) >30 kms^"*^, compared to 
17% of the Control sample stars.) However, a Kolmogorov-Smirnov (K-S) test comparing the two 
samples suggests that the difference is not significant; the test gives a probability of 0.74, indicat- 
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ing that the v sin{i) distributions of the Periodic and Control samples do not differ. Therefore the 
idea that studies of pre-main-sequence rotation are biased by the use of photometrically-measured 
rotation periods is not borne out by our data. 

T Tauri stars of similar mass and age arc observed to have a wide range of rotation periods 
and, in some cases, a bimodal period distribution (Stassun ct al. 1999; Herbst et al. 2000). A 
proposed explanation for the bimodality involves the concept of "disk-locking" , in which a magnetic 
interaction between the disks around T Tauri stars and the stars themselves regulates their rotation 
rates (Konigl 1991; Shu et al. 1994; Ostriker & Shu 1995). In this picture, kilogauss-strength 
magnetic fields intercept the circumstellar disk, channel accretion flow onto the star's surface, and 
"lock" the star to the disk. The rotation rate of the star is then forced to be the same as that of 
the disk at the locking radius. The specifics of disk-locking, e.g., how long it lasts, why some stars 
may be locked to their disks while others are not, and exactly what role it plays in a star's angular 
momentum evolution, axe unclear. Moreover, the observational evidence for disk-locking occurring 
during the T Tauri phase has recently been under debate (Stassun et al. 1999; Herbst et al. 2000). 

In the disk-locking picture, a star that is rotating relatively slowly is thought to be locked 
to its circumstellar disk (or only recently unlocked), whereas a fast rotator either was never disk- 
locked or has had time to spin up since it became unlocked. As a test of this idea, another useful 
comparison to make would be to compare the v sin{i) distributions of stars with disks and without 
them. To construct a sample of stars likely to have disks, we chose objects from our sample with 
characteristics like those of Classical T Tauri Stars (CTTS). CTTS often have strong emission lines 
and are thought to be experiencing accretion onto their surfaces from a circumstellar disk, resulting 
in surface hot spots that cause irregular light variations and flaring (Bertout 1989). Weak T Tauri 
Stars (WTTS), on the other hand, often have low-amplitude (<0.5 mag), periodic light variations, 
possibly due to the presence of large cool spots on their surfaces, and lack obvious signatures of an 
accretion disk (Herbst et al. 1994). 

We used the observational diagnostics from Hillenbrand et al. (1998) to choose a sample of 
possible CTTS and WTTS. These diagnostics, namely the near-infrared excess A(/c — K) and 
Ca II emission line strength W\{Ca II), are determined by Hillenbrand et al. and applied to stars 
in the ONC to investigate the frequency of circumstellar disks in the cluster. Here we adopt a 
combination of the criteria used by Hillenbrand et al. The disked (CTTS) sample we chose consists 
of stars that have no measured period, and have either (1) A(/c — K) > 0.3 and Ca II emission, 
or (2) A(/c — K) > 0.5. The non-disked (WTTS) sample is comprised of stars with measured 
periods, and either (1) A{Ic-K) < 0.1 and Wx{Ca II) > 0, or (2) A{Ic-K) < 0. We tried to use 
fairly stringent criteria while still choosing a good-sized sample of objects for each category, but 
the WTTS and CTTS samples are nevertheless small in size, with 30 and 36 objects, respectively. 
The V sin{i) distributions for the WTTS and CTTS samples are shown in Figure 10. Again, the 
numbers in each bin are normalized by the number of stars in the sample. The distributions appear 
to differ in the predicted sense; that is, the CTTS sample (stars likely to have disks) does contain 
a slightly larger fraction of slow rotators compared to the WTTS sample. A K-S test comparing 
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the two distributions, however, yields a probability of 0.26, meaning that any difference between 
them is not statistically significant. 

3.1.2. Comparing v sin{i) to Velocity Estimated from Rotation Period 

Because we know the rotation periods of most of the stars for which we measured vsin{i), 
we can compare the equatorial rotation velocities (calculated using the periods and radii estimated 
from the luminosities and effective temperatures) for those stars with their measured v sin{i) values. 
This comparison is shown in Figure 11: v sin{i) values for 153 stars are plotted versus equatorial 
velocities derived from the equation Veq = 2TrR/P. (Note that one star — JW 834 — has a v sin{i) 
and period, but has neither a luminosity nor an effective temperature determination, so its radius 
cannot be calculated.) SB2 candidates arc shown with open circles. For stars with v sin{i) less than 
or equal to our estimated limit of 11 kms^^, upper limit symbols are plotted. The solid line in the 
figure marks v sin{i) = v^g, and the dotted line marks v sin{i) = {TT/4)veq. Assuming randomly- 
oriented stellar axes, the observed mean value of sin{i) should be 7r/4, or 0.785 (Chandrasekhar & 
Miinch 1950). Therefore, at least in theory, the points are expected to scatter about the dotted 
line. 

Figure 11 shows first of all that there is a strong correlation between v sin{i) and Veq, thus 
demonstrating conclusively that the periodicity of T Tauri stars is caused by the rotation of stars 
with spotted surfaces. The correlation of these quantities is strong: results from two different 
non-parametric tests indicate that the probability that there is no correlation is between 10^^^ 
and 10~^^. In addition, the v sin{i) = v^q limit is respected by nearly all the stars to within 
the errors of the measurement. The one notable exception, JW 275, is discussed below. This is 
a gratifying result and an important confirmation of the basic assumptions underlying pre-main- 
sequence rotation studies to date. 

A noticeable feature of Figure 11 is that the majority of points fall below the v sin{i) = 
{'!r/4)veq line, indicating that the average sin{i) value we observe is lower than the expected value. 
To investigate the significance of this deviation, we begin by calculating a mean sin{i) including 
all 153 stars with both v sin{i) and Veq values. This yields (sin(i)) = 0.82 it 0.09, with a standard 
deviation of 1.1. This result is dominated by a few very high values of sin{i) that cannot simply 
be understood in terms of random errors. The most extreme example is JW 275, which has a 
radius of 2.3 Rq and WO period of 20.1 days, implying an equatorial velocity of 5.9 kms~^. 
Its measured vsin{i), however, is 78 kms^^, completely inconsistent with expectation. This star 
represents a clear example of a gross error; either its rotation period or its v sin{i) is incorrect. The 
star's period was measured from only a single longitude, and therefore could be a beat period. If 
the measured 20.1-day period is actually the beat period between a 1-day sampling interval and 
a 1.05-day rotation period, then Veq for the star should be 112 kms~^. Combining this with the 
measured v sin{i) yields the quite reasonable sin{i) value of 0.69. Rotation periods close to one 
day are difficult to distinguish from beat periods in data obtained at a single longitude. Clearly 
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this star and others like it need to be removed from the sample before drawing inferences about 
the observed distribution of sin{i). Stars whose v sin{i) measurements are too small compared to 
their likely errors should also be removed before {sin{i)) is calculated. 

Accordingly, we first restrict the sample to those stars with v sin{i) greater than 11 kms~^, our 
approximate measurement limit based on the discussion in Section 2.4. Wc also remove from the 
sample the remaining stars (which include JW 275) with predicted v^q values less than 11 kms~^. 
The rationale for eliminating those stars is first of all that if the period and radius are accurate and 
Veq is less than ~11 kms""^, we do not expect to measure a reliable v sin{i) given our resolution. 
Alternatively, if the period and radius are in error, the star should for that reason be culled from 
the sample. For the latter reason we exclude three stars with sin{i) > 1 that may be cases similar 
to JW 275. Applying these conditions leaves a sample of 86 stars with precise sin{i) estimates. 
This sample yields (sin(i)) = 0.64ib0.02 with standard deviation = 0.23. (If we include the above- 
mentioned three stars with sin{i) > 1, the result is the same within the errors: (sin(i)) = 0.67±0.03, 
standard deviation 0.26.) Furthermore, if we construct an even more robust sample of stars with 
V sin{i) and v^q greater than 13.6 kms~^, our estimate of the formal spectral resolution, then 
{sin(i)) = 0.59±0.025 (standard deviation = 0.21). The average inclination for this sample of 68 
objects is 38±2 degrees.^ 

Results similar to ours have been reported in previous studies that combined v sin{i) data with 
rotation periods of pre-main-sequence stars. Hartmann ct al. (1986) compared v sin{i) and V(,q for 
six T Tauri stars and found that the v sin(i) values of two of the slower rotators in their sample fell 
below the expected Veq values. Weaver (1987) combined data from Hartmann et al. (1986) with 
that from other studies to produce a sample of 20 stars for which both v sin{i) and period were 
known. Weaver found that the mean inclination for the sample of 19 stars with sin{i) < 1 was 40 
degrees rather than the expected value of 57 degrees. The average sin{i) value from Weaver's data 
is 0.63. More recently, Sodcrblom et al. (1999) observed 35 stars in the open cluster NGC 2264. 
For the ~20 stars in their sample with measured v sin{i) and period and sin{i) < 1, the average 
sin{i) value was 0.6, also in agreement with our result. 

There are a number of possible explanations for the discrepancy between observed average 
sin{i) values and theoretical expectation. To clarify the possible reasons, we express sin{i) in 
terms of quantities either observed directly or derived from observations by writing: 

sin{i) = [constant) [v sin{i)) P Teff^ L~^l'^ (1) 

Systematic errors in vsin{i), luminosity, effective temperature, and/or rotation period could give 
rise to the observed deviation from {sin{i)). The approximate size of the systematic error(s) would 



*We note that the errors on the mean values given here were calculated assuming a Gaussian distribution. This is 
not strictly valid as sin{i) is not expected to be normally-distributed. Our actual observed sin{i) values, however, do 
follow a smooth, fairly symmetric distribution with a prominent peak. Thus we have used the mean and its standard 
error to characterize the observed distribution. 
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be 0.2/0.6, or about 30%. Alternatively, the smaller-than-expected {sin{i)) might be due to a real 
physical phenomenon. We have explored each of these possibilities and address them below. 

Physical causes. — There may be an actual physical cause behind the lower average sin{i) we 
observe. For example, Weaver (1987) suggested that it could be the result of selection against 
stars with high axial inclination in rotation studies, because circumstellar material obscures the 
stellar photosphere. We have explored this idea and find from numerical simulations that in order 
to reduce the average value of sin{i) from 0.79 to 0.60, all stars with inclinations between 52 and 
90 degrees must be systematically excluded from period studies. Given a random distribution of 
rotation axes, this means about 60% of the stars would have to be excluded, which seems unlikely. 
Alternatively, it is possible that rotation axes in young clusters like the ONC and NGC 2264 are 
aligned rather than oriented randomly. Some have suggested that the formation of a star cluster 
from a rotating molecular cloud threaded by a uniform magnetic field should eventually lead to 
alignment of stellar rotation axes (see Bodenheimer 1995 and references therein). Evidence with 
regard to such alignment is mixed, with some data showing alignment of structures within molecular 
clouds (see Heiles et al. 1993 and included references) and some indicating that stellar rotation axes 
are randomly oriented in certain clusters (e.g., the Pleiades; Stauffer 1991). Another explanation 
could be that the prominent surface spots that give rise to T Tauri stars' photometric variability are 
more likely to appear near the poles of the stellar rotation axis than near the equator. Consequently 
the stars that are detected in photometric period studies (and thus appear in Figure 11) would tend 
to have more pole-on inclinations, and thus have sin{i) values that are truly lower than expected 
from a random distribution. On the other hand, if a star were seen exactly pole-on, there would 
be no photometric variation observed. In any case, Figure 9 shows that there is not a significant 
bias between the v sin[i) distributions of stars selected from period studies (the Periodic sample) 
compared to those selected via only their luminosities and effective temperatures (the Control 
sample). 

Systematic errors in vsin{i). — One way to produce a smaller-than-expected average sin{i) 
value is if our measured v sin{i) values are systematically low. Increasing our values by at least 
30% would bring {sin{i)) close to its expected value. The comparison of our v sin{i) values for 
the 14 stars we have in common with other studies (see Figure 8) shows, however, that for over 
half the stars, our measured sin[i) values are actually slightly larger than the values from other 
work. The overall mean difference between our v sin[i) values and those of the two other studies is 
0.5±2.0 kms^^. This, coupled with the fact that other independent studies of pre-main-sequence 
stars in NGC 2264, Taurus and Orion also found average sin{i) values of approximately 0.6, makes it 
seem unlikely that a systematic v sin{i) error is wholly responsible for the {sin{i)) discrepancy. On 
the other hand, it is difficult to rule out the possibility that some smaller systematic error in v sin{i) 
is contributing to the problem. In Section 3.1 we explored systematic uncertainties associated with 
line blending and possible mismatches in spectral type between the targets and cross-correlation 
templates. In those cases the associated uncertainties were accounted for by our estimated errors. 
But there could be other systematic effects at work. For example, narrow-lined spectra from dwarf 
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stars are used as cross-correlation templates against which the pre-main-sequence stellar spectra 
are compared. It is feasible that a difference in the limb darkening law between dwarfs and pre- 
main-sequcncc stars could have an effect on the v sin{i) calibration, introducing a systematic error. 
Such a systematic error could potentially affect the other v sin{i) studies mentioned, most of which 
also used dwarfs as templates for the cross-correlation. 

Errors in the photometric periods. — Still another explanation for the lower average sin{i) 
value could be that the rotation periods inferred from photometric studies are, in some cases, not 
the actual rotation periods of the stars. Instead, they could be false detections, harmonics of the 
rotation period, beat periods between an observing interval (normally one day) and the rotation 
period, or, in the case of spectroscopic binaries, related instead to the orbital period of the system. 
We consider each of these ideas in turn. 

It is unlikely that many of the photometric periods are false detections rather than true periods, 
since statistical tests arc used to eliminate these and, in many individual cases of low sin[i) values, 
there are confirming period determinations at different epochs (see Herbst ct al. 2000). Moreover, 
the presence of false detections should produce a random scatter in the v sin{i)—Veq plane, tending 
to obscure the observed correlation rather than skew it towards small values of sin{i). Harmonics, 
on the other hand, could explain at least some of what is displayed in Figure 11. T Tauri stars 
such as V410 Tau (Herbst 1989) and some stars in the ONC (Stassun et al. 1999; Herbst et al. 
2001) have displayed so-called "period doubling", i.e., the period determined at one epoch is later 
shown to be half of the true rotation period because the star originally had significant spots on two 
opposing hemispheres. The effect of correcting a period-doubling error would be to decrease Vgq by a 
factor of two for any such star. Taking the 68-star sample described above and arbitrarily doubling 
the period of every star with an inferred sin{i) < 0.60 (34 of the 68 stars) increases {sin{i)) to the 
expected value. Whether this is part or all of the problem can only be determined by continued 
careful monitoring of ONC stars to see whether period doubling occurs for a significant number 
of them. Comparisons made so far indicate that this is not a major contributor to the problem 
(Herbst etal. 2001; Carpenter, Hillenbrand & Skrutskie 2001). 

Another explanation for those stars with larger-than-expected Vgg could be that the observed 
photometric periods are actually beat periods. This would preferentially affect stars with shorter 
periods (i.e., ^1 day or so) since the periods are closer to the observing frequency. Both the 
WO (Herbst et al. 2000) and ESO (Herbst et al. 2001) data sets were obtained at single longi- 
tudes (although with multiple observations per night) and are thus more susceptible to the beat 
phenomenon than the Stassun et al. (1999) study, which combined data from widely-separated 
longitudes. Herbst et al. (2001) compared their data with Stassun et al. (1999) and found that 
of 111 stars monitored by both studies, 98 yielded identical periods, two were likely examples of 
period doubling, one appeared to be a false detection, and the other 10 were related as beat peri- 
ods with a 1-day sampling interval. This suggests that about 10% of the sample of periodic stars 
might be contaminated by beat periods or harmonics. Again, the continued monitoring of ONC 
stars, preferably from observatories well separated in longitude, is necessary in order to determine 
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whether these phenomena are important in the interpretation of Figure 11. 

Lastly, it is possible that some of the periodicity detected in photometric studies is of orbital 
origin. Close spectroscopic binaries are variable for a number of reasons including eclipses, tidal 
distortion and the reflection effect. Such stars may not be rotationally locked and observed periods 
could be the orbital period or a related harmonic. Alternatively, in the case of tidal locking, both 
tidal distortion and reflection can lead to period doubling. As mentioned in Section 3.1, seven 
stars show evidence indicating that they might be spectroscopic binaries. Five of these stars have 
measured periods and are marked with open circles in Figure 11. Of the five stars with the largest 
Veq values in the figure, two are SB2 candidates. If these stars' observed periods arise from orbital 
motion rather than rotation, then they should not be plotted on Figure 11. 

In summary, despite the fact that the rotation periods are more precise than the other quan- 
tities contributing to sin{i) (typically being known to ~1% or better; sec, e.g., Choi & Herbst 
1996), they are prone to gross errors such as aliasing and period doubling. These effects probably 
account for the small number of outliers removed from the sample used to calculate our observed 
{sin{i)). It is difficult to see how they could account for the remaining systematic effect, which is 
accompanied by a relatively small scatter in sin{i). In order for such period errors to account for 
the remaining systematic effect, a significant number of rotation periods would need to have been 
underestimated. This could be the case if period-doubling plays a major role. This seems unlikely, 
however, in light of the general agreement found between periods measured by different observers 
at different epochs with different sampling intervals (e.g., Stassun et al. 1999, Herbst et al. 2001, 
Rebuh 2001, Carpenter etal. 2001). 

Overestimated luminosities. — Equation (1) shows that sin{i) depends on L~^/^; as a result, 
systematically overestimated luminosities for our target stars could give rise to a lower-than- 
expected {sin{i)). For the sample of 68 stars, the true stellar luminosities would have to be 
approximately 60% of the current measured values in order to yield an average sin{i) consistent 
with the expected value. Weaver (1987) has suggested that the observed luminosities of prc-main- 
scquence stars might be artificially inflated due to the presence of a circumstellar disk. This idea 
could explain why the sin{i) discrepancy apparent in our data has been seen by others (namely 
Soderblom et al. 1999 and Weaver 1987), who have observed pre-main-sequence stars in different 
clusters with different observational methods. On the other hand, this requires that the large ma- 
jority of our target stars have a relatively luminous disk, i.e., one that is responsible for ~40% of 
the total optical luminosity of the disk-star system. Hillenbrand etal. (1998) estimate the disk 
fraction in the ONC to be between 55% and 90%, so the majority of ONC stars may indeed have 
disks. However when we examine the A(Ic — K) values of the 68-star sample, we find that only 
about 40% have A(Ic — K) > 0.3, indicating the presence of a prominent disk. There are uncer- 
tainties associated with determining A (7c — K), and Hillenbrand et al. note that the errors (which 
are expected to be primarily negative, reducing A(/c — K) from its true value) may be as large 
as 0.1—0.3 magnitudes for stars with spectral types K2 to M3. Nevertheless it seems unlikely that 
the large majority of our stars have disks that are sufficiently luminous to account for the observed 



-20- 



{sin{i)) discrepancy. 

Another uncertainty associated with the luminosities is the extinction correction used to derive 
them. Systematically overestimated extinction values would result in too-large luminosities. Our 
luminosities come from Hillenbrand (1997), who suggests that if anything the luminosity values 
may be underestimated, if there is unaccounted-for accretion activity or a nonstandard or spatially- 
variable extinction law in the ONC. In calculating the overall uncertainty on the luminosities, 
Hillenbrand assumes an error in Ay of 0.5 magnitudes, which is larger than expected and which 
corresponds to an uncertainty of 0.12 in log(L*/L0). Other sources of error that may contribute 
are intrinsic uncertainties due to photometric variability, errors from translating spectral type and 
color to luminosity, and uncertainties in the stellar distances. Hillenbrand estimates that these 
contribute uncertainties of 0.12, 0.15 and 0.075 dex, respectively, for a total combined uncertainty 
of less than 0.2 dex, on average. Subtracting this amount from the luminosities for the sample of 
68 stars yields {sin{i)) = 0.75±0.03, which is close to the expected value. But again, Hillenbrand 
suggests that if a systematic error is present in the luminosities, it would likely go the other way, 
leading to underestimated values instead of too-large ones. 

A detailed discussion of the observational uncertainties associated with estimates of stellar 
luminosities is given in Hartmann (2001), who explores the effect of such errors on age estimates 
in star- forming regions. In addition to discussing each of the sources of error mentioned above, 
Hartmann suggests that unresolved binary companions can have a potentially large effect on lu- 
minosity determinations. Moreover, this would be a systematic effect: the presence of unidentified 
companion stars around our targets would bias their luminosities to higher values. Hartmann mod- 
els the effect of unresolved binaries on luminosities and estimates that the shift in log L would 
be ~0.2, under the assumption that nearly all stars have unidentified companions. As explained 
above, a shift of this magnitude could bring {sin{i)) close to its expectation value, so this idea 
could potentially help explain the observed discrepancy. 

Underestimated effective temperatures. Because sin(i) oc TJj, the sin{i) values we derive are 
much more sensitive to errors in the effective temperatures than to errors in the luminosities. The 
effective temperatures we use come from Hillenbrand (1997), who determined spectral types for 
ONC stars and then converted them to T^s values using the calibration of Cohen &; Kuhi (1979), 
with a few modifications. We find that increasing T^s for our sample stars by 400 to 600 Kelvin 
(without changing the luminosity) yields good agreement between the measured and expected sin{i) 
distributions. Figure 12 shows the effect of a 600-K increase in temperature in the v sin{i) versus 
Veq plane. All 153 stars that have v sin{i) and v^q measures arc shown. The points now scatter 
around the dotted line representing v sin{i) = {'K/4:)veq- The figure demonstrates that increasing 
the temperatures can rectify the {sin{i)) discrepancy; however, this fix of the problem requires that 
the temperatures have been significantly underestimated. 

An examination of the literature suggests that it is possible that uncertainty inherent in the 
conversion from spectral type to T^e could cause the temperatures to be off by the required amount. 
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Padgett (1996) and King (1998), for example, discuss problems with the existing spectral type— Tgff 
calibration for pre-main-sequence stars. These include the uncertainty with regard to whether a 
dwarf or giant temperature scale is more appropriate, as discussed by Hillenbrand (1997) and 
White et al. (1999). King estimates that the size of the errors in T^s determinations for pre-main- 
sequence stars is 500—700 K. Padgett finds differences as large as 960 K between T^s values derived 
from line ratio measurements in high resolution spectra and those obtained from spectral type- 
Teg calibration. She finds that temperatures derived from metallic line ratios are generally higher 
than those determined from lower-resolution spectral classification. The assumption that pre-main- 
sequence stellar atmospheres can be precisely calibrated by comparison with the atmospheres of 
dwarfs or evolved stars is likely to break down at some point. A major difference is that pre- 
main-sequence stars have spots — indicating the presence of strong magnetic fields — covering a 
large portion of their surfaces. The effect of these spots on the spectral type-Teg calibration is not 
well-understood. Moreover, Cohen & Kuhi (1979) note that it is unclear exactly what the physical 
meaning of Tpff is in a star with different emitting regions (such as the photosphere, chromosphere, 
extended atmosphere, and circumstellar disk) that each contribute significant amounts of flux. 

To summarize, we flnd a significant departure of the distribution of measured sin{i) values from 
what is expected for a random distribution of observed inclination angles. This departure has also 
been seen in previous studies of pre-main-sequence stars. It could have a number of explanations, 
including real physical conditions, selection effects, and/or systematic errors in one or more of the 
quantities that go into calculating sin{i). If this were entirely due to an error in the spectral type-to- 
Teff calibration, an adjustment by 400 to 600 Kelvin towards hotter temperatures than determined 
by Hillenbrand (1997) is required to solve the problem. This is large, but perhaps not impossible, 
given the current state of our understanding of effective temperatures of pre-main-sequence stars. If 
the Tefj values of the stars have been underestimated, the masses and ages (determined by plotting 
the stars on an H-R diagram and interpolating the masses and ages from evolutionary model tracks) 
also need to be adjusted. (For example, a 0.25-Mq star becomes a O.4-M0 star if its temperature is 
increased by 600 K.) This in turn would have important implications for conclusions drawn about 
pre-main-sequence stellar evolution in the ONC and other clusters and associations. 



3.2. Stellar Radii 

A main goal of this study is to combine the v sin{i) measurements we obtain with rotation 
period measurements from other studies in order to constrain the radii of the stars in our sample. 
Because P v sin{i) = 2 ir R sin{i), the v sin{i) and period of a star can be used to calculate its 
minimum radius, R sin{i). With enough data, the rotational velocities and periods can provide 
statistical estimates of radii (without the sin{i) ambiguity) for stars grouped by their similar 
characteristics. Direct measurements of the radii of pre-main-sequence stars are hard to come by, 
since so far only four pre-main-sequence eclipsing binaries have been discovered (Popper 1987; 
Casey et al. 1995, 1998; Andersen 1991; Mamajek et al. 2000; Covino et al. 2000). Therefore 
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other methods for independently estimating stellar radii can be valuable tests of the reliability of 
the radii of pre-main-sequence stars and hence the accuracy with which they can be placed on the 
theoretical H-R diagram. 

We have combined period and v sin{i) measurements to calculate minimum radii for individual 
stars with v sin{i) greater than our estimated limit of 11 kms~^. It should be emphasized that 
these minimum radii are not subject to the uncertainties in luminosity and effective temperature 
discussed in the previous section, but instead rely only on the measured period and v sin{i). There 
are 101 such stars not including JW 275, which has a period that is likely in error, as noted in 
Section 3.1.2. (Calculating a minimum radius using JW 275's measured period of 20.1 days yields 
an R sin{i) of 30.8 R©; if its true period is instead 1.05 days, its R sin{i) is 1.61 R©.) Minimum 
radii for the 101 stars are given in Table 1 in units of the solar radius, and Figure 13 shows a 
histogram of the values. They range from less than one to more than five solar radii, with the 
majority (80%) having R sin{i) < 2.5 R©. The median of the distribution is 1.7 R©. 

Next we turn to making statistical estimates of the radii of groups of stars located in the same 
part of the H-R diagram. This calculation requires a sample of stars with reliable v sin{i) and 
period measurements, as well as known L and T^s values. L and T^g are needed only so that the 
stars can be grouped in a way that will make the radii derivation meaningful; again, the radii we 
calculate are completely independent of these two quantities. The same potential sources of error 
that were relevant when we explored the magnitude of the {sin{i)) discrepancy in Section 3.1.2 
must be considered here as well. For example, stars with sin(i) values much larger than 1 may 
have rotation periods that are in error. In addition, v sin{i) values near our estimated measurement 
limit (11 kms~^) may be less certain than larger values. To derive the most robust average radii 
possible with these data, we use the 68-star sample — consisting of stars with accurate periods and 
V sin{i) >13.6 kms^^ — that was constructed to help characterize the {sin{i)) problem. 

We want to calculate average radii for groups of stars that are in approximately the same 
evolutionary stage. Consequently the stars must be binned according to their location in the 
theoretical H-R diagram. Even if there is a systematic error in T^s such as that discussed in 
Section 3.1.2, the relative locations of the sample stars should be fairly constant, so that grouping 
them by effective temperature should still be reasonable. Using as a guide the way in which 
the stars are arranged on the H-R diagram, a series of equal-sized bins was chosen that would 
maximize the number of stars within each bin (and the number of stars that could be binned) 
while still retaining adequate resolution to derive meaningful average radii. The bins are shown as 
solid lines in Figure 14. They each cover 0.048 units in log Tpff space and 0.37 units in log (L/Lq) 
space. The 68 sample stars used here are shown as filled circles; three of them fall outside of the 
chosen bins. The dotted lines marked on the figure are lines of constant radius, at intervals of 1 R©, 
beginning at the bottom with 1 R©. 

Eleven of the 30 bins contain at least two stars. For these bins, an average stellar radius 
was calculated, by first computing R sin{i) for each individual star using the measured v sin{i) 
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and period, averaging the projected radii for the stars in the bin, and finally dividing the average 
R sin{i) by the average expected value of sin{i), 7r/4. The radii are given in Table 3. The first 
three columns of the table specify the bin number and luminosity and effective temperature range 
spanned by the bin. Column 4 gives the number of stars in the bin. Column 5 gives the mean 
R sin{i) value and its error, and column 6 lists the mean radius and error. The last column lists, 
for comparison, the average radius corresponding to the center of each L and Tes bin. The mean 
difference between the radii calculated from v sin{i) and P and the radii for the centers of the bins 
is —0.70, with a standard deviation of 0.74. Some difference between the R values calculated using 
V sin{i) and period and those based on the locations of the stars in the H-R diagram is of course 
expected, in light of the data shown in the v sin{i) versus Vgq plot and the accompanying discussion. 

The calculation of {R) from period and v sin{i) yields a measurement of the radii of our sample 
stars that is independent of the usual methods for deriving stellar radii, and provides a snapshot of 
the evolution of low-mass stars as they contract toward the main sequence. To investigate whether 
we can actually observe that stars are contracting (i.e., use our data to show this geometrically), 
we have grouped the bins in Table 3 according to their T^q values, then ordered bins with the same 
Tcff values by decreasing luminosity. Figure 2 shows that the evolutionary tracks in the portion 
of the H-R diagram in which our stars are located run fairly vertically. Therefore, examining the 
average R sin(i) values in sets of vertical bins (i.e., comparing bin Al to A2, Bl to B2, etc.) should 
tell us whether stars of similar mass are actually contracting as they approach the Zero Age Main 
Sequence (ZAMS). 

Examination of Table 3 shows that the mean radii do decrease with decreasing luminosity 
for bins in the same temperature range. However, the relative errors on [R sin{i)) are as large as 
20—30% in some cases. As a result, there is evidence for contraction at greater than the three-sigma 
level in only one pair of bins, B2 to B4. This is the only instance in which non-adjacent bins can be 
compared to each other. So although the mean values do decrease, we cannot state with absolute 
certainty that we are observing pre-main-sequence contraction for stars within a given temperature 
range. If we could include more stars with reliable v sin{i) and period in each bin (and by doing 
so reduce the errors on [R sin(i))) or could compare bins that span a larger range in radius, we 
might find more definitive evidence for contraction of stars in a specific temperature range. 

We can look for overall evidence of contraction for the entire sample of 68 stars by combining 
the information in the five sets of vertical bins shown in Table 3. For each set of vertical bins (A, B, 
C, D, and E) listed in the table, we can calculate the relative decrease in {R sin{i)) with decreasing 
luminosity. For example, the relative change in radius (A R sin{i)/R sin{i)) for bins Al to A2 
is 0.21. Calculating this quantity for each of the five sets of vertical bins and averaging the result 
yields a mean relative change in projected radius of 0.49it0.16. This result is positive — i.e., shows 
that stars in our sample are contracting as they approach the main sequence — at the three-sigma 
level. This is the first time that pre-main-sequence contraction has been verified geometrically. 

To conclude, these data show marginally significant evidence supporting a real spread in radius 
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(and therefore age) at a given mass in the ONC. It is clear that to proceed further (for example, 
determining an accurate age spread for the cluster) requires a larger sample of stars with accurate 
period and v sin(i) measurements. However, the promise of this method has now been clearly 
demonstrated. 

4. Summary & Conclusions 

We have obtained high-dispersion spectroscopic data for 256 pre-main-sequence stars in the 
Orion Nebula Cluster, with the aim of measuring their v sin{i) values. Our intention was to 
combine the data with that from previous studies of the ONC that measured rotation periods via 
photometric monitoring, in order to better understand the rotational evolution of T Tauri stars 
and to arrive at statistical estimates of the sizes of these stars. A set of 118 objects with known 
periods (the Periodic sample) was chosen along with 138 stars from the same portion of the H-R 
diagram whose periods were not known (the Control sample). Most of the targets are low- mass 
ONC members with spectral types from early-K to mid-M. Spectra with sufficient signal-to-noise 
to measure v sin{i) were obtained for 238 of the 256 targets. Our results are as follows: 

1. Approximately two-thirds of the 238 sample stars with good-quality spectra have v sin{i) less 
than 20 kms~^, and the v sin{i) values of the remaining stars are distributed in a tail extending to 
~100 kms-^ 

2. We have compared the v sin{i) distributions of the Periodic and Control samples in order to 
determine whether photometric rotation period studies are truly representative of the rotational 
properties of pre-main-sequence stars. The differences between the distributions are not statistically 
significant. The same is true when we compare a sample of stars likely to be accreting material 
from their circumstellar disks (CTTS), and perhaps "disk-locked", to a sample that does not show 
evidence for accretion (WTTS). 

3. We find a strong correlation between v sin{i) and equatorial rotation velocity v^q derived by 
combining the stars' periods with radii calculated from their published L and Teff values. This 
shows definitively that the observed periodicity of T Tauri stars is, in the majority of cases, caused 
by the rotation of spots on their surfaces. 

4. The average sin{i) value for our stars, calculated from sin{i) = v sin{i)/veq, is significantly 
lower than the mean value expected for a random distribution of stellar rotation axes. This result, 
which has been seen in previous studies with smaller samples, could be due to systematic errors in 
one or more of the observed quantities that go into the sin{i) calculation (namely v,sin{i), period, 
L, and Tcff), or to a real physical phenomenon. We explore these possibilities and find that the 
correct sin{i) value is produced if we assume that the effective temperatures of our program stars 
have been underestimated by 400—600 Kelvin. 

5. We have calculated minimum radii for stars that have rotation period measurements and v sin{i) 
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greater than our estimated measurement limit. The R sin{i) values range from 0.4 to 5.3 Rq, with 
four-fifths of the sample falling between 0.5 and 2.5 Rq. 

6. We have derived average radii for groups of stars with similar masses (i.e., in similar locations on 
the H-R diagram). The mean radii for stars within the same mass/temperature range do decrease 
as the stars get closer to the ZAMS, but in most cases the uncertainties in the mean radii are too 
large to show definitively that we have observed contraction within a given mass range. However, we 
find that the mean relative change in radius with luminosity for the five different mass/temperature 
ranges is 0.49ib0.16, meaning that we do find evidence for pre-main-sequence contraction in the 
overall sample of stars used at the three-sigma level. 

In conclusion, combining v sin{i) measurements with photometric periods can yield valuable 
insight into the astrophysics of T Tauri stars. The methods used here allow us to test the lumi- 
nosities and effective temperatures of T Tauri stars and thereby test their masses and ages, since 
these are usually derived by comparing L and Tpg with pre-main-sequence evohitionary models. 
Furthermore, even with a modest number of stars with good-quality v sin(i) values and periods, we 
have shown geometrically for the first time that the stars in our sample are undergoing contraction 
as they evolve toward the main sequence. Applying the method to even larger samples of stars 
with well-determined v sm(z)'s and periods is likely to improve our understanding of this important 
stage of low-mass stellar evolution. 
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Fig. 1. — The locations on the sky of the 256 ONC stars observed spectroscopically for this study 
are shown here as filled circles. 1053 stars included in the Jones & Walker (1988) proper motion 
survey are shown as small dots. The four Trapezium stars are located at the center of the small 
box in the middle of the field. 
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Fig. 2. — The positions on the H-R diagram of the 256 stars included in this study are shown here. 
Filled circles are stars in the Periodic sample, and open circles are Control sample stars. Luminosity 
and effective temperature values are from Hillenbrand (1997). Pre- main-sequence evolutionary 
tracks are from D'Antona &; Mazzitelli (1994, model 1). 
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Fig. 3. — Examples of template star and target star spectra obtained for this study. The top 
spectrum is from Gliese 144, a K2V star used as a spectral template in the cross-correlation pro- 
cedure. Its rotational velocity is below our v sin{i) measurement limit. The bottom four panels 
show spectra of Jones & Walker (1988) stars with varying projected rotational velocities; the stellar 
absorption lines become broader and more shallow with increasing vsin(i). Emission lines from 
the Orion Nebula appear in some of the spectra. 
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Fig. 4. — Calibration curve for the narrow-lined template star Gliese 114, created by cross- 
correlating the original spectrum of this star against a series of artificially-broadened versions of 
the spectrum. The FWHM of the cross-correlation peak was measured and is plotted here against 
the velocity corresponding to each broadened spectrum (points) . The polynomial fit to the FWHM 
vs. velocity data is shown as a dotted line. The relationship between the width of the peak and 
V sin{i) is fairly linear above 7 kms~^. 
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Fig. 5. — Shown here are the peaks of the cross-correlation functions for the target stars JW 790 
(top panel) and JW 669 (bottom two panels). The cross-correlation functions are produced by 
cross-correlating the stars against the narrow-lined template star closest in effective temperature 
to the target. For JW 790, the cross-correlation peak is regular in shape and well-approximated 
by a Gaussian function. The peaks for JW 669, made using spectra taken in January 1997 and 
December 1997, appear more like two overlapping Gaussians, an indication that this star may be 
a spectroscopic binary. 
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Fig. 6. — Plotted here are the cross-correlation function peaks for six of the seven double-lined 
spectroscopic binary candidates found in our data set. (The other SB2 candidate, JW 669, is 

shown in Figure 5.) These cross-correlation peaks show structure (sometimes appearing as multiple 
overlapping Gaussians) and are wider than typical peaks, possibly due to the presence of two sets 
of spectral lines in the target star spectra. 
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Fig. 7. — Histogram of the heliocentric radial velocity values for the stars for which v sin{i) was 
measured. The mean and standard deviation of the distribution are 26.7 kms~^ and 5.6 kms~^, 
respectively. 




Fig. 8. — Comparison of v sin{i) values for 14 stars observed for this study and values from Duncan 
(1993) (squares) and Wolff, Strom, & Hillenbrand (2001) (triangles). 
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Fig. 9. — The top plot shows the distribution of v sin{i) values for the 238 stars for which v sin{i) 
was successfully measured. The bottom two plots are the v sin{i) distributions for the Periodic and 
Control sample stars. Although there appears to be a slight excess of rapid rotators in the Periodic 
sample, a K-S test comparing the v sin{i) distributions of the Periodic and Control samples shows 
that they are not significantly different from each other. 
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Fig. 10. — Comparison of the v sin{i) distributions for a sample of 36 stars likely to have accretion 
disks (CTTS) and 30 stars not likely to have disks (WTTS). The fraction of slow rotators is slightly 
higher for the CTTS sample, but a K-S test shows that the distributions do not differ significantly. 
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rotational velocity (2ttR/P) [km/s] 



Fig. 11. — Measured v sin[i) versus equatorial rotation velocity for the 153 stars for which the 
latter quantity can be calculated. SB2 candidates are plotted with open circles, and v sin{i) values 
less than or equal to 11.0 kms~^ are shown as upper limits. Equatorial velocities are derived from 
the equation Veg = 2ttR/P. The solid line marks v sin{i) = Veq, and the dotted line marks v sin{i) 
= {Tr/4)veq, where (7r/4) is the expected average value of sin{i) for a random distribution of stellar 
rotation axes. 
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rotational velocity (SttR/P) [km/s] 



Fig. 12. — A new version of Figure 11: measured v sin{i) versus equatorial rotation velocity is 
shown for 153 stars. Here, the equatorial velocities have been calculated using Teg values that are 
600 K hotter than published values. The points scatter around the dotted v sin{i) = {'ir/4)veq line, 
as expected for a group of stars with randomly-oriented rotation axes. 
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Fig. 13. — Minimum radii for 101 stars with measured period and v sin{i) > 11.0 km s ^, calculated 
from the equation R sin{i) = P v sin{i) / (277) . 80% of the stars have R sin{i) < 2.5 R0. 
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Fig. 14. — This figure shows the locations of the radii bins described in Section 3.2. The bins are 
marked with soUd fines; each is 0.048 units wide in log Tgs and 0.37 units wide in log (L/Lq). The 
bin numbers are designated in the figure; e.g., the bin in the bottom right corner is bin E6. 68 
stars with v sin{i) > 13.6 kms^^ and known rotation period, luminosity and efi^ective temperature 
appear as filled circles. Dotted lines mark lines of constant radii in the log-L, log-Tcfj plane, ranging 
from 1—10 R0 in increments of 1.0. For bins containing two or more stars, average radii have been 
calculated and are given in Table 3. 



Table 1. ONC Stars Observed for this Study 



JW 


Sam 


Mem 


I 


V-I 


log(L) 


log(Teff) 


A(I-K) 


TyA(Call) 


Per 


Src 


V sin{i) 


R sin{i) 


SB2? 


Date Obs 






(%) 


(mag) 


(mag) 


(L©) 


(K) 


(mag) 


(A) 


(days) 




(km s-l) 


(R©) 






2 


C 


97 


12.2 


1.61 


0.26 


3.643 


0.42 








<11.0 






Jan 97,Dec 97 


3 


C 





13.0 


2.50 


0.24 


3.580 


0.55 


2.0 


3.43 


E 


29.8 ± 4.5 


2.02 




Dec 97 


15 


C 


99 


14.1 


2.66 


-0.18 


3.562 


0.11 


1.7 


9.56 


E 


<11.0 






Dec 97 


17 


c 


99 


13.5 


3.00 


-0.20 


3.509 


0.10 


0.8 


3.14 


E 


19.9 ± 3.9 


1.23 




Dec 98 


18 


c 


99 


14.4 


2.66 


-0.46 


3.544 


0.54 


0.6 






13.1 ± 3.3 






Dec 97 


20 


c 


99 


13.8 


2.67 


-0.40 


3.518 


0.10 


0.7 


0.67 


E 


70.3 ± 11.6 


0.93 




Dec 98 


22 


c 


99 


13.9 


2.38 


-0.34 


3.555 


0.14 


2.2 






<11.0 






Dec 97 


25 


c 


99 


15.2 


2.54 


-0.88 


3.500 


0.52 


0.0 


2.28 


E 


18.0 ± 4.0 


0.81 




Dec 98 


27 


c 


93 


11.1 


1.31 


0.72 


3.695 


0.22 








38.7 ± 2.2 






Dec 97 


29 


c 


97 


11.7 


1.60 


0.67 


3.695 


0.32 








23.8 ± 1.4 






Jan 97 


31 


c 


99 


14.2 


2.44 


-0.58 


3.526 


0.01 


1.4 






<11.0 






Dec 98 


36 


c 


99 


14.3 


3.08 


-0.30 


3.526 


0.46 


-4.1 


2.72 


E 


21.9 ± 3.9 


1.18 




Dec 98 


37 


c 


92 


12.9 


1.68 


0.14 


3.679 


0.62 


-8.4 


8.27 


E 


14.9 ± 0.9 


2.43 




Jan 97 


39 


c 


99 


14.3 


2.40 


-0.63 


3.526 


0.14 


0.8 






11.7 ± 1.1 






Dec 98 


46 


c 


97 


11.1 


1.27 


0.73 


3.708 


0.10 








14.6 ± 0.6 






Dec 97 


50 


c 


98 


11.7 


1.75 


0.44 


3.602 


1.37 


-14.6 






16.4 ± 2.3 




Y 


Dec 98 


51 


c 


99 


14.4 


2.37 


-0.57 


3.555 


1.05 


0.0 






12.4 ± 1.5 






Dec 98 


55 


c 


95 


14.4 


2.67 


-0.60 


3.526 


0.36 


1.0 






<11.0 






Dec 97 


63 


c 


99 


12.8 


2.06 


0.26 


3.623 


0.20 


1.9 






17.4 ± 1.2 






Jan 97 


65 


p 


99 


13.8 


2.25 


-0.24 


3.580 


0.50 


0.0 


7.39 


V 


12.0 ± 1.1 


1.75 




Dec 97 


73 


c 


99 


14.0 


1.85 


-0.59 


3.562 


1.65 


-34.0 












Dec 98 


75 


p 


99 


11.1 


1.61 


0.89 


3.695 


0.09 




3.45 


V 


29.6 ± 1.8 


2.02 




Jan 97 


76 


p 


99 


13.6 


2.97 


-0.23 


3.509 


0.16 


1.3 


6.33 


V 


<11.0 






Jan 97 


77 


p 


99 


13.6 


2.16 


-0.21 


3.580 


0.42 


1.8 


1.50 


EVS 


38.8 ± 5.3 


1.15 




Dec 97 


83 


p 


99 


14.7 


3.12 


-0.29 


3.544 


0.37 


-9.2 


7.72 


V 


18.3 ± 5.9 


2.79 




Dec 97 


91 


c 


99 


13.5 


2.44 


-0.24 


3.509 


-0.06 


1.2 


17.08 


E 


12.9 ± 1.2 


4.35 




Dec 97 


95 


c 


95 


12.8 


1.77 


-0.12 


3.591 


0.22 


0.0 






61.5 ± 9.0 






Dec 97 


98 


c 


99 


14.6 


3.20 


-0.51 


3.509 


0.17 


0.0 


2.34 


E 


20.9 ± 2.8 


0.97 




Dec 98 


99 


p 


99 


13.1 


2.51 


0.13 


3.562 


0.50 


1.4 


1.70 


V 


60.2 ± 10.5 


2.02 


Y 


Jan 97,Dcc 97 


106 


c 


99 


13.4 


2.42 


-0.29 


3.526 


0.10 


1.5 


1.70 


E 


23.0 ± 2.9 


0.77 




Dec 98 


111 


c 


98 


13.0 


1.83 


-0.15 


3.591 


0.24 


2.2 


4.94 


E 


<11.0 






Jan 97 


114 


p 


99 


15.3 


2.14 


-1.01 


3.562 


2.01 


-1.4 


8.73 


V 


13.8 ± 4.6 


2.38 




Dec 97 


116 


p 


99 


12.1 


1.44 


0.20 


3.643 


0.09 


1.6 


2.34 


V 


39.0 ± 3.3 


1.80 




Jan 97 


117 


c 


99 


13.2 


2.13 


-0.05 


3.580 


0.77 


-6.3 


8.87 


E 


<11.0 






Jan 97 


122 


c 


98 


14.4 


3.02 


-0.54 


3.494 


-0.17 


0.0 


0.98 


E 


62.2 ± 17.1 


1.20 


Y 


Dec 97 


124 


c 


99 


13.9 


2.76 


-0.43 


3.518 


-0.02 


0.8 






14.7 ± 3.8 






Dec 97 


126 


c 


99 


14.4 


2.55 


-0.69 


3.526 


0.01 


0.8 


2.78 


E 


18.4 ± 2.8 


1.01 




Dec 97 


127 


c 


99 


14.3 


2.94 


-0.45 


3.518 


1.05 


-2.8 






13.9 ± 2.1 






Dec 97 


128 


p 


97 


12.6 


2.82 


0.28 


3.535 


0.32 


0.0 


8.83 


V 


<11.0 






Dec 97 


133 


p 


99 


13.2 


2.08 


-0.25 


3.555 


0.26 


1.6 


2.02 


V 


30.9 ± 4.5 


1.23 




Dec 97 


135 


p 


99 


14.4 


2.93 


-0.46 


3.526 


1.11 


-2.3 


3.69 


V 


36.8 ± 8.5 


2.68 




Dec 97 


136 


p 


99 


13.8 


2.81 


0.06 


3.571 






8.65 


V 


11.6 ± 1.7 


1.98 




Dec 98 


140 


c 


99 


13.9 


3.10 


-0.19 


3.518 


-0.29 


3.8 


4.58 


E 


<11.0 






Dec 98 


146 


c 


99 


13.5 


2.27 


-0.33 


3.544 


0.32 


0.0 






11.2 ± 1.3 






Dec 98 


148 


p 


99 


15.4 


3.73 


0.07 


3.580 


0.08 


1.6 


1.37 


VS 








Dec 97 


149 


p 


99 


15.6 


3.25 


-0.70 


3.526 


0.24 


-23.7 


2.83 


V 








Dec 97 


150 


c 


99 


14.0 


2.87 


0.27 


3.623 


1.13 


-1.7 






<11.0 






Dec 98 


152 


c 


99 


14.1 


2.78 


-0.42 


3.526 


0.39 


2.6 












Dec 98 
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Table 1 — Continued 



JW 


Sam 


Mem 


I 


V-I 


log(L) 


log(Tpff) 


A{I-K) 


Wx (Call) 


Per 


Src 


V sin(i) 


R sin(i) 


SB2? Date Obs 






(%) 


(mag) 


(mag) 




(K) 


(mag) 


(A) 


fdavsl 




(km s-i) 


(Rn) 




157 


P 


99 


10.2 


1.56 


1.19 


3.679 


0.11 




17.40 


V 


<11.0 




Dec 98 


158 


P 


99 


13.3 


2.15 


-0.23 


3.555 


0.71 


1.5 


1.95 


V 


19.5 ± 2.8 


0.75 


Jan 97 


159 


P 


99 


13.9 


2.46 


-0.37 


3.544 


-1.72 


0.3 


1.12 


V 


76.0 ± 18.6 


1.68 


Dec 97 


163 


C 


99 


11.8 


1.30 


0.22 


3.643 


0.53 


1.6 






<11.0 




Dec 97 


164 


C 


99 


14.0 


2.38 


-0.47 


3.544 


0.63 


-0.7 


6.52 


E 


12.2 ± 1.6 


1.57 


Dec 98 


165 


P 


97 


11.7 


1.82 


1.49 


3.897 


0.11 




5.77 


V 


17.7 ± 1.1 


2.02 


Dec 98 


166 


C 


99 


14.2 


2.81 


-0.46 


3.494 


-0.26 


0.0 


0.67 


E 


29.6 ± 8.3 


0.39 


Dec 97 


167 


P 


99 


14.9 


3.08 


-0.03 


3.602 


0.86 


1.4 


3.43 


E 


<11.0 




Jan 97 


169 


P 


99 


15.9 




-1.30 


3.526 






3.85 


V 






Dec 98 


174 


P 


99 


14.1 


2.40 


-0.59 


3.535 


0.12 


3.7 


1.36 


EVS 


46.2 ± 5.4 


1.24 


Jan 97,Dec 


175 


P 


99 


14.1 


2.41 


-0.55 


3.526 


0.56 


0.0 


9.19 


V 


<11.0 




Dec 97 


176 


C 


99 


11.6 


1.89 


0.47 


3.591 


0.54 


0.8 






<11.0 




Dec 98 


180 


c 


99 


14.4 


2.94 


-0.54 


3.494 


0.05 


0.2 






41.9 ± 9.6 




Dec 98 


187 


c 


99 


14.3 


2.85 


0.22 


3.643 


0.53 


1.4 


14.41 


E 


<11.0 




Dec 98 


191 


p 


99 


14.2 


2.78 


-0.37 


3.535 


0.18 


0.0 


8.63 


V 


<11.0 




Dec 97 


203 


c 


98 


13.8 


2.48 


-0.11 


3.580 


0.63 


1.0 






<11.0 




Dec 97 


208 


c 


96 


14.3 


2.52 


-0.56 


3.535 


0.24 


0.0 






<11.0 




Dec 98 


211 


c 


99 


14.0 


2.61 


-0.42 


3.535 


0.24 


1.4 


5.46 


E 


13.2 ± 2.3 


1.42 


Dec 97 


220 


p 


99 


12.7 


2.25 


0.07 


3.555 


0.67 


1.0 


2.33 


V 


19.7 ± 3.1 


0.91 


Jan 97 


221 


c 


98 


11.3 


2.66 


1.41 


3.679 


-0.01 


1.5 






56.9 ± 5.3 




Jan 97 


222 


p 


98 


14.1 


2.85 


-0.23 


3.544 


0.71 


1.2 


5.17 


V 


14.8 ± 3.0 


1.51 


Dec 97 


225 


c 


99 


13.4 


2.30 


-0.18 


3.555 


0.14 


0.5 






25.8 ± 3.5 




Jan 97 


229 


c 


99 


14.2 


2.41 


-0.43 


3.555 


0.37 


2.1 


7.75 


E 


12.9 ± 1.8 


1.98 


Dec 97 


232 


p 


79 


11.7 


1.86 


0.82 


3.695 


0.23 




5.08 


V 


34.0 ± 2.2 


3.41 


Dec 97 


237 


c 


99 


15.3 


3.22 


-0.78 


3.509 


0.75 


0.0 


5.22 


E 


<11.0 




Dec 98 


239 


p 


99 


12.7 


2.25 


0.05 


3.555 


0.38 


1.5 


4.45 


V 


23.5 ± 3.5 


2.07 


Jan 97 


243 


p 


99 


13.6 


2.36 


-0.31 


3.544 


0.64 


1.2 


10.20 


V 


12.5 ± 1.8 


2.52 


Dec 97 


248 


p 


99 


12.7 


2.17 


0.10 


3.571 


0.59 


-0.9 


6.86 


V 


12.7 ± 1.6 


1.72 


Jan 97 


249 


c 


99 


13.6 


2.35 


0.35 


3.695 


1.14 


1.7 






25.9 ± 1.7 




Dec 97 


250 


p 


99 


14.0 


3.00 


-0.24 


3.526 


0.27 


0.9 


2.71 


V 


20.5 ± 2.7 


1.10 


Dec 97 


254 


p 


99 


14.3 


3.21 


-0.24 


3.526 


0.50 


5.6 


3.52 


V 


19.4 ± 3.2 


1.35 


Dec 98 


256 


c 


99 


14.1 


1.13 


-0.52 


3.518 




2.1 






<11.0 




Dec 98 


272 


c 


99 


14.1 


2.30 


-0.46 


3.555 


0.72 


-2.2 






<11.0 




Dec 98 


275 


p 


99 


13.0 


2.28 


-0.13 


3.544 


0.76 


1.8 


20.10 


V 


77.6 ± 15.7 




Dec 98 


284 


p 


99 


14.6 


3.11 


-0.41 


3.526 


1.18 


-2.0 


3.11 


EV 






Jan 97 


291 


c 


99 


12.9 


2.07 


-0.06 


3.562 


0.32 


1.5 


1.91 


E 


14.3 ± 1.9 


0.54 


Dec 98 


305 


c 


99 


14.6 


2.39 


-0.72 


3.518 


0.42 


0.0 






13.8 ± 2.5 




Dec 98 


308 


p 





12.6 


2.54 


0.06 


3.526 




0.7 


4.17 


V 


30.2 ± 4.8 


2.49 


Dec 98 


310 


c 


98 


13.8 


3.57 


-0.05 


3.500 


0.58 


0.8 






54.5 ± 15.2 




Dec 98 


311 


p 


99 


14.9 


3.13 


-0.60 


3.518 


0.22 


0.0 


6.14 


V 


<11.0 




Dec 98 


315 


p 


99 


14.8 


2.85 


-0.73 


3.518 


-0.03 


0.0 


8.85 


V 


<11.0 




Dec 97 


320 


p 


99 


14.7 


2.37 


-0.08 


3.695 


0.98 


0.3 


5.41 


V 


11.6 ± 2.0 


1.24 


Dec 97 


321 


c 


95 


13.5 


2.09 


-0.19 


3.580 




-5.1 


7.44 


E 


13.8 ± 1.3 


2.03 


Dec 97 


328 


c 


99 


12.9 


1.64 


0.09 


3.661 


1.02 


0.0 


18.83 


E 


14.3 ± 1.2 


5.32 


Dec 97 


330 


p 


99 


11.8 


1.42 


0.45 


3.679 


0.08 


1.8 


1.57 


V 


25.5 ± 1.3 


0.79 


Jan 97 


332 


c 


99 


14.2 


2.07 


-0.59 


3.518 


2.46 


-0.9 






<11.0 




Dec 98 


333 


c 


99 


14.7 


1.92 


-0.56 


3.483 


0.54 


3.7 






<11.0 




Dec 98 


334 


p 


99 


13.6 


2.35 


-0.20 


3.562 


1.20 


-4.8 


5.27 


V 


16.9 ± 3.2 


1.76 


Dec 97 
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JW 


Sam 


Mem 


I 


V-I 


log(L) 


log(Tpff) 


A{I-K) 


Wx (Call) 


Per 


Src 


V sin{i) 


R sin(i) 


SB2? Date Obs 






(%) 


(mag) 


(mag) 




(K) 


(mag) 


(A) 


fdavsl 




(km s-i) 


(Rn) 




335 


C 


99 


15.6 


1.31 


-1.21 


3.562 


3.10 


-26.0 










Dec 98 


337 


P 


99 


11.7 


1.76 


0.34 


3.591 


0.52 


1.9 


19.50 


V 


<11.0 




Dec 97 


345 


P 


81 


12.4 


2.38 


0.37 


3.571 


0.42 


0.0 


8.21 


V 


32.5 ± 4.9 


5.27 


Dec 97 


347 


P 


99 


12.5 


1.75 


0.09 


3.602 


0.81 


0.0 


7.33 


vs 


<11.0 




Dec 97 


348 


c 


99 


11.2 


1.53 


0.92 


3.719 


1.00 








21.7 ± 1.1 




Dec 98 


352 


p 


99 


10.4 


2.10 


1.09 


3.591 




1.9 


8.00 


V 


14.3 ± 0.8 


2.26 


Jan 97 


362 


p 


99 


13.9 


3.23 


-0.29 


3.500 


0.38 


0.0 


2.73 


V 


23.2 ± 8.3 


1.25 


Dec 97 


363 


c 


98 


12.7 


2.51 


0.12 


3.544 


0.40 


1.1 






33.4 ± 7.6 




Jan 97 


365 


p 


99 


11.7 


1.96 


0.82 


3.679 


0.40 




4.08 


V 


41.6 ± 3.0 


3.35 


Jan 97 


373 


c 


99 


13.5 


3.37 


0.99 


3.679 


0.20 


2.1 


9.81 


E 


16.3 ± 0.8 


3.16 


Dec 98 


378 


p 


99 


12.9 


2.11 


0.05 


3.580 


1.08 


1.5 


9.03 


V 


<11.0 




Dec 97 


379 


p 


99 


15.2 


3.87 


-0.58 


3.483 


0.07 




11.30 


VS 


<11.0 




Dec 97 


381 


p 


99 


13.2 


2.11 


-0.14 


3.571 


1.26 


0.0 


16.20 


V 


<11.0 




Dec 97 


388 


p 


99 


13.6 


2.62 


0.22 


3.602 


0.06 


1.2 


9.08 


V 


11.9 ± 0.9 


2.13 


Dec 97 


398 


c 


99 


14.5 


2.44 


-0.41 


3.580 


0.30 


2.2 






<11.0 




Dec 98 


406 


p 


99 


13.9 


3.41 


0.37 


3.562 


0.31 


1.7 


2.25 


V 


30.1 ± 6.6 


1.34 


Dec 97 


417 


p 


98 


13.9 


2.44 


-0.22 


3.571 


0.56 


0.0 


7.41 


V 


<11.0 




Dec 97 


421 


c 


99 


11.5 


1.46 


0.47 


3.643 


0.80 


-3.5 






<11.0 




Dec 98 


429 


c 


99 


12.7 


1.75 


0.57 


3.753 




4.5 






13.2 ± 0.7 




Dec 98 


433 


c 


93 


11.9 


1.59 


0.43 


3.661 


0.21 








<11.0 




.. Jan 97,Dec 


437 


p 


99 


11.8 


2.07 


0.56 


3.602 


0.53 


1.4 


2.34 


V 


45.7 ± 4.3 


2.11 


Jan 97 


439 


p 


99 


15.0 


3.33 


-0.12 


3.562 


0.62 


0.7 


8.30 


V 


<11.0 




Dec 98 


440 


c 


99 


14.5 


3.09 


-0.44 


3.518 


-0.35 


0.0 






49.3 ± 8.3 




Y Dec 98 


443 


c 


99 


14.0 


0.95 


-0.55 


3.544 


0.84 


-1.5 






11.5 ± 2.2 




Dec 98 


454 


p 


99 


11.9 


2.10 


0.78 


3.661 


0.39 


-1.0 


8.67 


V 


20.5 ± 1.2 


3.51 


Dec 97 


460 


c 


95 


11.7 


1.58 


0.65 


3.695 


0.88 








40.9 ± 4.2 




Dec 98 


466 


p 


99 


14.4 


2.88 


-0.56 


3.518 


0.35 




6.17 


V 


12.7 ± 3.0 


1.55 


Dec 97 


470 


p 


99 


12.5 


2.63 


0.84 


3.661 




6.1 


10.70 


V 


23.9 ± 2.0 


5.05 


Dec 98 


478 


p 


99 


11.6 


2.28 


0.69 


3.580 


0.43 




5.74 


V 


35.3 ± 3.5 


4.00 


Dec 97 


479 


p 


99 


11.0 


1.80 


1.16 


3.719 


0.57 


2.0 


8.71 


V 


25.4 ± 1.4 


4.37 


Dec 98 


481 


p 


99 


14.7 


2.99 


-0.52 


3.526 


0.35 


0.0 


6.56 


V 






Dec 97 


482 


c 


99 


13.3 


2.38 


0.10 


3.591 


0.60 


-50.3 






13.9 ± 3.0 




Dec 98 


485 


p 


99 


14.3 


2.49 


-0.42 


3.555 


1.18 


0.8 


9.66 


V 


<11.0 




Dec 97 


487 


c 


99 


12.8 


1.78 


0.07 


3.623 


0.28 


1.8 






12.5 ± 0.9 




Jan 97 


498 


p 


99 


13.8 


2.08 


-0.44 


3.562 






7.42 


V 


11.7 ± 1.4 


1.72 


Dec 97 


501 


c 


99 


13.2 


1.82 


-0.28 


3.580 


0.62 


1.1 






<11.0 




Dec 98 


502 


c 


99 


13.5 


2.57 


-0.22 


3.535 


0.97 


0.0 






23.8 ± 4.7 




Dec 98 


526 


p 


98 


11.9 


1.90 


0.52 


3.623 


0.74 


1.2 


1.68 


V 


95.3 ± 14.6 


3.16 


Jan 97 


527 


c 


99 


13.0 


2.09 


-0.06 


3.571 


0.63 


0.8 






12.7 ± 1.5 




Dec 98 


534 


c 


85 


12.7 


1.43 


-0.03 


3.544 


1.33 


-69.5 










Dec 98 


544 


p 


99 


10.8 


2.22 


1.30 


3.670 


-0.17 


1.6 


3.42 


V 


58.8 ± 5.8 


3.97 


Jan 97 


565 


c 


99 


14.5 


2.95 


-0.62 


3.509 


0.22 


0.0 






<11.0 




Dec 98 


567 


p 


99 


9.7 


1.73 


1.43 


3.661 


-0.17 


-1.7 


8.53 


V 


22.2 ± 1.1 


3.74 


Dec 97 


576 


p 


99 


13.1 


2.28 


-0.07 


3.555 


0.23 


0.0 


2.01 


V 


30.5 ± 4.5 


1.21 


Dec 97 


585 


c 


99 


11.5 


1.37 


0.39 


3.602 




1.8 






18.4 ± 0.9 




Dec 98 


589 


p 


99 


11.3 


2.12 


1.23 


3.719 


0.11 


-0.6 


14.30 


V 


14.6 ± 0.7 


4.12 


Dec 98 


592 


c 


99 


14.4 


2.58 


-0.57 


3.535 


0.06 


0.0 


7.92 


E 


<11.0 




Dec 98 


601 


p 


98 


12.5 


1.25 


0.05 


3.679 


0.05 


2.5 


2.22 


V 


11.9 ± 0.5 


0.52 


Jan 97 
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JW 


Sam 


Mem 


I 


V-I 


log(L) 


log(Tpff) 


A{I-K) 


Wx (Call) 


Per 


Src 


V sin{i) 


R sin(i) 


SB2? Date Obs 






(%) 


(mag) 


(mag) 




(K) 


(mag) 


(A) 


fdavsl 




(km s-i) 


(Rn) 




625 


C 


99 


UA 


1.07 


-0.73 


3.562 


1.37 


6.1 










Dec 98 


629 


C 


99 


13.5 


4.34 


1.47 


3.643 


1.41 


-16.9 






56.3 ± 11.3 




Dec 98 


631 


c 


99 


12.9 


2.20 


0.22 


3.602 


1.29 


1.6 


0.96 


E 


39.7 ± 3.7 


0.75 


Dec 98 


636 


p 


99 


13.4 


2.46 


0.04 


3.575 


0.79 




5.59 


V 


<11.0 




Jan 97 


639 


p 


99 


13.6 


2.53 


-0.37 


3.526 






5.24 


V 


12.7 ± 1.7 


1.31 


Dec 97 


641 


p 


99 


11.5 


1.74 


0.91 


3.719 


0.12 


0.0 


3.17 


V 


46.8 it 3.8 


2.93 


Jan 97 


647 


c 


99 


13.7 


2.31 


-0.22 


3.494 




-5.0 






13.8 ± 1.3 




Dec 98 


648 


p 


99 


13.2 


2.58 


0.24 


3.580 




0.9 


34.50 


V 


<11.0 




Dec 98 


661 


c 


99 


13.5 


1.57 


-0.39 


3.580 


1.80 


2.1 






<11.0 




Dec 97 


662 


c 


95 


14.4 


1.78 


-0.57 


3.623 


1.36 


-1.0 






30.7 ± 4.8 




Dec 98 


664 


p 


99 


14.1 


2.88 


0.26 


3.623 


-0.18 


1.4 


7.20 


V 


<11.0 




Dec 97 


669 


p 


99 


10.8 


1.67 


1.12 


3.708 


0.14 


1.8 


1.81 


V 


58.4 ± 4.0 


2.09 


Y Jan 97,Dec 


673 


p 


99 


13.1 


1.83 


0.06 


3.643 


0.81 


1.3 


1.44 


EV 


18.7 ± 1.3 


0.53 


Jan 97 


678 


p 


99 


12.3 


1.84 


0.55 


3.695 


0.24 




13.00 


V 


12.4 ± 0.7 


3.18 


Dec 97 


683 


p 


99 


11.8 


1.85 


0.53 


3.623 


0.28 


1.9 


11.50 


V 


<11.0 




Dec 98 


694 


c 


99 


14.7 


2.99 


0.05 


3.623 


0.18 


-3.8 






<11.0 




Dec 98 


701 


c 


99 


12.4 


1.26 


0.25 


3.722 


1.27 


-8.3 






12.4 ± 0.8 




Dec 98 


704 


p 


99 


14.7 


2.41 


-0.59 


3.562 






6.87 


V 


12.7 ± 3.7 


1.72 


Dec 97 


706 


c 


99 


11.8 


1.34 


0.39 


3.679 


-0.03 








12.0 ± 0.6 




Dec 98 


710 


p 


99 


12.7 


1.17 


-0.04 


3.562 


0.68 


-2.1 


7.59 


V 


<11.0 




Dec 98 


716 


p 


99 


14.4 


2.88 


-0.57 


3.500 


1.15 


-2.7 


3.95 


V 


<11.0 




Dec 98 


717 


p 


99 


16.1 




-1.39 


3.535 




0.0 


4.00 


V 






Dec 98 


719 


c 


99 


14.5 


3.06 


-0.56 


3.494 


0.10 


0.0 


6.88 


E 


<11.0 




Dec 97 


721 


p 


99 


13.1 


2.08 


-0.06 


3.580 


0.82 


-3.8 


2.45 


V 


<11.0 




Dec 97 


727 


p 


99 


13.9 


2.44 


-0.38 


3.544 


0.75 


0.8 


6.03 


V 


13.3 ± 2.7 


1.58 


Dec 97 


728 


c 


99 


13.6 


2.32 


-0.20 


3.562 


0.54 


0.9 


13.87 


E 


<11.0 




Dec 97 


731 


p 


88 


12.3 


1.74 


0.10 


3.591 


0.28 


2.2 


7.64 


V 


12.7 ± 0.9 


1.92 


Dec 97 


733 


c 


99 


13.1 


2.33 


0.04 


3.571 


0.58 


1.7 


3.28 


E 


20.8 ± 3.0 


1.35 


Dec 98 


736 


c 


99 


13.9 


2.27 


-0.25 


3.580 


0.65 


1.6 


11.68 


E 


<11.0 




Dec 98 


744 


c 


99 


13.4 


1.75 


-0.36 


3.562 


2.42 


-0.9 






<11.0 




Dec 98 


751 


c 


99 


13.5 


2.47 


-0.07 


3.562 


0.78 


-3.0 






<11.0 




Dec 97 


756 


c 


99 


12.8 


2.15 


0.09 


3.580 


1.67 


-2.0 






33.2 ± 2.7 




Dec 98 


758 


p 


99 


14.1 


3.45 


-0.15 


3.509 


0.48 


-2.0 


2.51 


V 


16.0 ± 3.0 


0.79 


Dec 98 


765 


p 


97 


13.2 


2.67 


-0.18 


3.522 


0.09 


1.3 


2.42 


V 


14.2 ± 2.4 


0.68 


Dec 97 


779 


c 


95 


12.5 


1.34 


0.05 


3.661 


0.55 








18.7 ± 1.1 




Dec 98 


782 


c 


99 


14.2 


2.14 


-0.41 


3.591 


2.18 


-2.5 






<11.0 




Dec 98 


786 


p 


89 


13.9 


2.39 


-0.23 


3.571 


1.30 


0.0 


8.81 


V 


<11.0 




Dec 97 


787 


p 


99 


12.7 




-0.07 


3.571 




0.6 


9.09 


E 


<11.0 




Dec 98 


788 


p 


99 


14.5 


3.48 


-0.28 


3.509 


0.28 


0.0 


4.82 


EV 


<11.0 




Dec 97 


790 


p 


79 


10.9 


1.06 


0.73 


3.719 


0.12 


2.2 


2.74 


V 


39.8 ± 2.8 


2.15 


Jan 97 


792 


p 


99 


13.6 


2.64 


0.08 


3.580 


1.05 


1.4 


10.30 


V 


<11.0 




Dec 97 


794 


p 





10.9 


0.81 


0.62 


3.722 


0.03 




2.62 


V 


<11.0 




Dec 98 


795 


p 


99 


10.6 


1.19 


1.02 


3.744 


0.12 


1.7 


1.55 


V 


60.9 ± 5.4 


1.86 


Jan 97 


811 


p 


99 


13.4 


2.56 


0.13 


3.580 


1.25 


0.0 


11.10 


V 


<11.0 




Dec 98 


813 


p 


99 


12.4 


1.75 


0.05 


3.580 


0.56 


1.9 


2.85 


EV 


26.0 ± 2.5 


1.46 


Dec 97 


815 


p 


99 


13.8 


2.69 


-0.34 


3.526 


0.52 


0.0 


6.42 


V 


12.9 ± 1.9 


1.64 


Dec 97 


817 


p 


99 


13.8 


3.04 


-0.30 


3.500 


0.22 




2.23 


EV 


40.0 ± 6.5 


1.76 


Jan 97,Dec 


818 


c 


98 


12.1 


1.82 


0.20 


3.591 


0.26 








14.5 ± 0.9 




Dec 98 
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JW 


Sam 


Mem 


/ 


V-I 


log(L) 


log(Tefr) 


A{I-K) 


H^A(Call) 


Per 


Src 


V sin(i) 


R sin{i) 


SB2? Date Obs 






(%) 


(mag) 


(mag) 


(Ln) 


(K) 


(mag) 


(A) 


(days) 




(km s-i) 


(R«) 




819 


P 


99 


13.2 


2.85 


-0.06 


3.518 


0.12 


1.3 


1.21 


E 


27.7 ± 9.4 


0.66 


Dec 97 


823 


P 


99 


13.9 


3.28 


-0.31 


3.494 


0.25 


0.0 


9.00 


V 


<11.0 




Dec 98 


826 


P 


91 


14.0 


1.38 


-0.60 


3.602 


3.01 


-6.3 


9.57 


V 


<11.0 




Dec 97 


830 


P 


99 


13.0 


2.36 


0.01 


3.555 


0.63 


1.8 


9.23 


V 


<11.0 




Dec 98 


834 


C 


99 


13.0 


2.24 










7.01 


E 


13.4 ± 1.2 


2.44 


Dec 98 


835 


P 


99 


14.1 


2.52 


-0.54 


3.526 


0.24 


0.0 


8.85 


V 


<11.0 




Dec 97 


836 


P 


99 


13.4 


2.20 


-0.11 


3.580 


1.10 


0.9 


12.20 


V 


<11.0 




Dec 97 


837 


P 


99 


13.3 


2.21 


-0.06 


3.580 


0.62 


1.2 


10.60 


V 


<11.0 




Dec 97 


839 


P 


99 


13.5 


3.70 


0.70 


3.562 


-1.41 


1.2 


7.52 


V 


<11.0 




Dec 97 


843 


P 


98 


13.8 


2.45 


-0.33 


3.544 


0.02 


0.0 


0.84 


E 


89.5 ± 16.3 


1.48 


.. Jan 97,Dec 97 


850 


P 


99 


13.4 


2.14 


-0.11 


3.591 


0.73 


-4.1 


7.78 


V 


<11.0 




Jan 97 


855 


P 


99 


14.1 


2.34 


-0.51 


3.544 


0.79 


0.0 


7.08 


V 


11.8 ± 2.0 


1.65 


Dec 97 


862 


c 


99 


15.5 


3.80 


-0.77 


3.471 


0.24 


-4.5 


3.49 


E 






Dec 97 


863 


p 


99 


14.0 


2.62 


-0.50 


3.518 


0.71 


0.0 


7.91 


V 


<11.0 




Dec 97 


864 


c 


99 


14.4 


2.46 


-0.70 


3.526 


0.59 


0.7 


3.70 


E 






Dec 98 


866 


p 


98 


11.8 


2.13 


0.99 


3.708 


0.18 




6.76 


V 


27.4 ± 1.6 


3.66 


Jan 97 


868 


c 


99 


13.2 


1.96 


0.16 


3.661 


0.42 


1.4 






13.7 ± 0.8 




Dec 98 


872 


p 


99 


15.2 


3.24 


-0.83 


3.494 


0.60 


0.4 


4.19 


V 






Dec 97 


880 


p 


99 


14.2 


2.74 


-0.51 


3.509 


-0.10 


0.0 


1.58 


E 


19.4 ± 1.9 


0.61 


Dec 98 


883 


c 


99 


13.3 


2.62 


0.12 


3.562 






0.85 


E 


61.3 ± 15.3 


1.03 


Y Jan 97 


887 


c 


97 


10.3 


1.42 


1.39 


3.771 


0.05 








66.2 ± 6.3 




Dec 98 


890 


p 


99 


13.0 


2.11 


-0.16 


3.555 


0.46 


1.6 


1.10 


V 


61.6 ± 13.3 


1.34 


Jan 97 


892 


p 


99 


14.0 


2.22 


-0.13 


3.623 




1.0 


8.56 


V 


12.0 ± 0.9 


2.03 


Dec 97 


906 


c 


99 


14.2 


2.45 


-0.58 


3.526 


0.36 


1.6 


10.66 


E 


<11.0 




Dec 97 


907 


p 


99 


12.8 


1.60 


0.17 


3.679 


-0.20 


1.8 


1.65 


V 


37.3 ± 2.8 


1.22 


Jan 97 


911 


c 


99 


12.7 


1.62 


0.15 


3.661 


1.01 








<11.0 




Dec 98 


912 


c 


99 


13.9 


2.54 


-0.48 


3.526 


0.48 


1.2 






<11.0 




Dec 98 


914 


p 


99 


13.2 


2.09 


-0.23 


3.555 


0.71 


0.0 


7.71 


V 


<11.0 




Dec 97 


919 


c 


97 


14.4 


2.93 


-0.45 


3.526 


-0.22 


1.4 






<11.0 




Dec 98 


926 


c 


96 


12.8 


1.93 


0.11 


3.602 


0.09 


1.8 


5.54 


E 


12.3 ± 1.0 


1.35 


Dec 98 


927 


c 


99 


14.0 


2.29 


-0.43 


3.509 


0.01 


1.0 


2.02 


E 


22.4 ± 5.7 


0.89 


Dec 97 


930 


p 


99 


14.1 


2.56 


0.04 


3.623 


-0.71 


5.2 


2.88 


V 


15.0 ± 1.7 


0.85 


Dec 97 


933 


p 


98 


13.7 


2.48 


-0.31 


3.509 


0.05 


1.0 


5.98 


V 


<11.0 




Jan 97 


938 


c 


99 


12.6 


2.70 


0.99 


3.695 


0.19 


1.3 


5.15 


E 


21.7 ± 2.6 


2.21 


Dec 97 


940 


c 


95 


12.0 


1.63 


0.25 


3.535 


-0.13 








61.6 ± 14.9 




Dec 98 


945 


c 


99 


12.5 


2.40 


2.29 


4.111 


0.58 












Dec 98 


949 


c 


90 


12.3 


2.55 


0.47 


3.562 


0.55 


0.0 






67.3 d= 18.8 




Dec 97 


961 


c 


98 


12.3 


1.71 


0.19 


3.612 


0.18 


1.9 


1.43 


E 


66.2 ± 7.2 


1.87 


Y Dec 98 


962 


c 


99 


13.3 


1.97 


0.00 


3.623 


0.80 


1.8 


9.56 


E 


<11.0 




Dec 98 


969 


c 


99 


14.4 


2.54 


-0.54 


3.544 


0.58 


-2.7 






14.0 ± 1.8 




Dec 98 


971 


c 


83 


11.5 


1.51 


0.40 


3.602 


0.19 


2.4 






42.7 ± 3.6 




Dec 98 


972 


p 


99 


13.3 


2.02 


-0.28 


3.562 


0.38 


1.2 


8.07 


V 


11.1 ± 2.1 


1.77 


Dec 97 


988 


c 


99 


14.8 


3.12 


-0.34 


3.544 


0.81 


-4.2 






19.9 ± 6.1 




Dec 97 


990 


c 


99 


14.1 


2.88 


-0.42 


3.500 


-0.14 


0.6 


2.16 


E 


21.4 ± 3.6 


0.91 


Dec 98 


995 


c 


98 


13.0 


2.11 


-0.14 


3.555 


0.22 


0.9 






<11.0 




Jan 97 


998 


c 


99 


13.2 


2.32 


-0.24 


3.535 


0.11 








12.8 ± 1.6 




Dec 98 


1002 


c 


99 


14.6 


3.08 


-0.41 


3.471 


0.51 


-18.6 










Dec 98 


1004 


c 


99 


13.8 


2.10 


-0.40 


3.562 


0.31 


2.0 






<11.0 




Dec 98 



Table 1 — Continued 



JW Sam Mcin / V-I log(L) Iog(T,,(i) A(I-K) W'A(Call) Per 

('/-;.) (mug) (mug) (L.J (K) (mag) (A) (days) 



Src 



(km s ''■) 



B .swi(/) SB2? 
(« 



Date Obs 



1006 


C 


99 


14.7 


2.92 


-0.65 


3.518 


2.21 


1.2 






<11.0 




Dec 98 


1007 


C 


99 


14.5 


2.50 


-0.72 


3.526 


-0.00 


1.1 






11.1 ± 1.5 




Dec 98 


1008 


C 


92 


12.7 


1.71 


-0.08 


3.571 


0.03 


2.0 






<11.0 




Jan 97 


1013 


C 


99 


14.6 


2.30 


-0.74 


3.526 


0.64 


-26.5 






25.6 ± 5.7 




Dec 98 


1016 


c 


88 


13.0 


1.92 


0.16 


3.643 


0.38 








12.4 ± 1.0 




Dec 98 


1020 


c 


99 


12.7 


1.40 


-0.08 


3.623 


-0.01 


1.7 






<11.0 




,Ian 97 


1021 


c 


36 


13.0 


1.91 


0.11 


3.633 


0.34 


1.6 


7.75 


E 


15.0 ± 1.7 


2.30 


Dec 97 


1026 


c 


98 


13.8 


2.24 


-0.27 


3.571 


1.41 


-1.3 


4.21 


E 


15.3 ± 1.7 


1.28 


Jan 97,Dec 


1030 


c 


99 


13.7 


2.18 


-0.43 


3.544 


0.36 


1.5 






11.9 ± 1.4 




Dec 98 


1033 


c 


80 


14.7 


2.04 


0.49 


3.917 


-0.31 


2.0 






34.2 ± 5.7 




Dec 98 


1035 


c 


98 


14.1 


2.35 


-0.55 


3.526 


0.09 


1.6 


0.70 


E 


77.2 ± 14.7 


1.07 


Dec 98 


1040 


c 


99 


13.6 


2.66 


-0.31 


3.518 


0.20 


0.6 






11.7 ± 2.1 




Dec 97 


1042 


c 


99 


16.0 


2.70 


-1.08 


3.483 


0.55 


0.0 










Dec 97 


1044 


c 


87 


14.6 


3.30 


-0.54 


3.500 


0.23 


-1.5 


2.98 


E 






Dec 98 


1047 


c 


99 


14.2 


2.84 


-0.51 


3.509 


-0.03 


0.6 






16.3 ± 2.8 




Dec 97 


1049 


c 


95 


11.8 


1.57 


0.30 


3.612 


0.10 


2.2 






18.6 ± 1.2 




Jan 97 



-48- 



Table 2. Gliese Stars Used as Narrow-lined Spectral Templates 



Gliese# 


Spectral Type 


V 


RA (1950) 


Dec (1950) 






(mag) 


(h m s) 


(deg ' ") 


75 


KOV 


5.6 


01 44 06 


+63 36 24 


144 


K2V 


3.7 


03 30 34 


-09 37 36 


114 


K7V 


8.9 


02 47 49 


+15 30 36 


15A 


M1.5V 


8.1 


00 15 31 


+43 44 24 


-111 


M2V 


7.5 


11 00 37 


+36 18 18 



-49- 



Table 3. Mean Radii 



Bin 


log L/LQ range 


log TeflF range 


N 


{Rsin{i)) 

(R©) 


{Rsin{i)) 
0.79 

(Ro) 


R for bin center 

(R©) 


Al 


1.11 


- 1.48 


3.682 - 


3.730 


3 


3.53 ± 0.72 


4.49 ± 0.92 


5.75 


A2 


0.74 


- 1.11 


3.682 - 


3.730 


5 


2.85 ± 0.32 


3.62 ± 0.41 


3.75 


Bl 


1.11 


- 1.48 


3.634 - 


3.682 


2 


3.86 ± 0.12 


4.91 ± 0.15 


7.17 


B2 


0.74 


- 1.11 


3.634 - 


3.682 


3 


3.34 ± 0.10 


4.25 ± 0.13 


4.68 


B4 


0.00 


- 0.37 


3.634 - 


3.682 


3 


1.18 ± 0.37 


1.51 ± 0.47 


2.00 


C3 


0.37 


- 0.74 


3.586 - 


3.634 


2 


2.64 ± 0.53 


3.36 ± 0.67 


3.81 


C4 


0.00 


- 0.37 


3.586 - 


3.634 


3 


1.16 ± 0.36 


1.48 ± 0.45 


2.49 


D4 


0.00 


- 0.37 


3.538 - 


3.586 


8 


1.52 ± 0.16 


1.94 ± 0.21 


3.11 


D5 


-0.37 


- 0.00 


3.538 - 


3.586 


10 


1.23 ± 0.11 


1.56 ± 0.14 


2.03 


E5 


-0.37 


- 0.00 


3.490 - 


3.538 


10 


1.08 ± 0.11 


1.37 ± 0.14 


2.53 


E6 


-0.74 - 


-0.37 


3.490 - 


3.538 


10 


0.92 ± 0.08 


1.17 ± 0.10 


1.65 



